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ABSTRACT 

The Herschel Space Observatory opens the sky for observations in the far infrared at high spectral 
and spatial resolution. A particular class of molecules will be directly observable; light diatomic 
hydrides and their ions (CH, OH, SH, NH, CH+, 0H+, SH+, NH+). These simple constituents are 
important both for the chemical evolution of the region and as tracers of high-energy radiation. If 
outflows of a forming star erode cavities in the envelope, protostellar far UV (FUV; Q < E^ < 13.6 
eV) radiation may escape through such low-density regions. Depending on the shape of the cavity, 
the FUV radiation then irradiates the quiescent envelope in the walls along the outflow. The chemical 
composition in these outflow walls is altered by photoreactions and heating via FUV photons in a 
manner similar to photo dominated regions (PDRs). 

In this work, we study the effect of cavity shapes, outflow density, and of a disk with the two- 
dimensional chemical model of a high-mass young stellar object introduced in the second paper in 
this series. The model has been extended with a self-consistent calculation of the dust temperature 
and a multi-zone escape probability method for the calculation of the molecular excitation and the 
prediction of line fluxes. 

We flnd that the shape of the cavity is particularly important in the innermost part of the envelope, 
where the dust temperatures are high enough (> 100 K) for water ice to evaporate. If the cavity shape 
allows FUV radiation to penetrate this hot-core region, the abundance of FUV destroyed species (e.g. 
water) is decreased. On larger scales, the shape of the cavity is less important for the chemistry in 
the outflow wall. In particular, diatomic hydrides and their ions CH"*", OH"*" and NH"*" are enhanced 
by many orders of magnitude in the outflow walls due to the combination of high gas temperatures 
and rapid photodissociation of more saturated species. The enhancement of these diatomic hydrides 
is sufRcicnt for a detection using the HIFI and PACS instruments onboard Herschel. The effect of 
X-ray ionization on the chemistry is found to be small, due to the much larger luminosity in FUV 
bands compared to X-rays. 

Subject headings: stars: formation - stars: individual: AFGL 2591 - molecular processes - ISM: 
molecules - X-rays: ISM 



1. INTRODUCTION 

Observations of young stellar objects (YSOs) en- 
velopes reveal a rich variety of morphologies (e.g. 
lArce k Sarge^ I2006L iJOrgensen et~a l. 2007). Many of 
the observed features are associated with outflows: For 
example, strong high-J CO emission along a cavity, 
etched out by the outflow, is foun d in the class I object 
HH 46 ()van Kempen et al.ll2009"bl ). They explain the CO 
narrow emission with a warm and FUV heated region, 
with FUV photons from either the bow shock, internal 
jet working surfaces or the accretion disk boundary layer. 
In this scenario, FUV radiation may travel freely in the 
low-densi ty cavity withou t being absorbed or scattered. 
Similarly, iBruderer et al.l (|2009a) have observed molecu- 
lar tracers for warm, dense and FUV irradiated gas along 
the outflow of the high-mass young stellar object AFGL 
2591. Since the line shape does not show clear signs of 
shocks, they suggest the FUV radiation to be of proto- 
stellar origin. Indeed, the hot photosphere of a luminous 
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young O or B star emits mostly at FUV wavelengths and 
a cavity that allows protostellar r adiation to escape ha s 
also been observed in that source (jPreibisch et al.ll2003D . 
This scenario of irradiated outflow walls along a cavity, 
swept out by an outflow, has also been used to explain 
other observations. For example, J0rgcnscn (2004) found 
emission of FUV enhanced CN along an outflow of a low- 
mass class source. Van Kempen et al. (,2009a ) explain 
the water maser emission at 183 GHz in Ser SMM 1 
with the increased water abundance in the warm, FUV 
or shock heated outflow walls. Outflow walls can also 
be sites of anomal o us molecula r excitation. For example 
iHog er heijde et al.l ()1998l 119991 ) observe strong emission 
of HON and HCO+ along the outflow of L 1527 and Ser 
SMM 1. These molecules have a large dipole moment 
and are efRciently excited by collisions with electrons, 
mixed fro m th e ioniz ed outflow into the outflow wall. 
Rawlings et al.l ()2004f ) explain the strong HCO+ emis- 
sion by shock liberation of molecules from ice m antles 
followed by photoprocessing. iSpaans et al.l (|1995[ ) model 
the strong ^'^CO 7 = 6^-5 emission in narrow lines ob- 
served toward several low-mass sources (e.g. TMR-1) by 
FUV heated outflow walls. Bruderer et al. (2009b) (sec- 
ond paper of this series, hereafter paper II) explain the 
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observed high abundance of the FUV enhanced molecule 
CO"*" in AFGL 2591 by a concave cavity that allows di- 
rect irradiation onto the outflow walls. 

Diatomic hydrides and their ions (XH, XH+) are 
corner-stones of the chemical network. For example wa- 
ter at high temperature (> 250 K) is mainly formed by 
the reaction of OH + H2. Due to their chemical prop- 
erties, discussed in this work, they trace dense, warm 
and FUV or X-ray irradiated gas like in outflow walls. 
Such tracers are valuable to constrain the FUV or X- 
ray radiation of embedded sources. This is of partic- 
ular interest for the innermost part of low-mass YSOs 
with protoplanetary disks being i onized and h eated by 
FUV radiation and X-rays (e.g. lBerginl[2009[ ). Here, 
we discuss the use of light diatomic hydrides with the 
common elements (X=0, C, N, S) as tracers of prot- 
stellar high-energy radiation. The chemistry of other di- 
atomic hydrides with heavier (e.g. X=Fe) and/or less 
abundant elements (e.g. X=F) is not discussed here, as 
their chemistry is poorly known and their use as tracers 
currently limited. However, with the ongoing search for 
new diatomic hydrides using the new Hcrschel Space Ob- 
servatory (e.g.Cernic haro et al. 2010), th is sitation may 
change. For example iNeufeld et al.l (|2010( ) report the de- 
tection of HF using the HIFI instrument onboard Her- 
schel. Due to the dissociation energy of HF being larger 
than of H2 , HF can be formed in an exothermic reaction 
of F with H2 and HF may become the dominant reservoir 
of fluorine (|Neufeld fc Wolfirell2009i) . The sihcon bear- 
ing SiH and SiH+ may be interesting in the context of 
shocked regions wh ere "sputtering" from dust grains (e.g. 
iFlower et al.lll996l ) can increase the gas phase abundance 
of silicon. 

The Herschel Space Observatory and Band 10 receivers 
of the upcoming Atacama Large Millimeter and sub- 
millimeter Array (ALMA) allow the study of important 
FUV and X-ray tracers, like diatomic hydrides, at high 
sensitivity and with good spatial and spectral resolu- 
tion for the first time. This raises the questions: What 
is the expected abundance and line flux of FUV en- 
hanced species in the scenario of directly irradiated out- 
flow walls? From which region of the envelope does the 
molecular radiation emerge and how are the molecules 
excited? How dependent are observable quantities (line 
fluxes or line ratios) on characteristics that cannot be 
constrained directly from observations (e.g. geometry in 
the innermost part of the envelope and chemical age)? 
To attack such questions, detailed numerical models of 
the chemistry of YSOs envelopes can be used. 

Models of the chemistry of YSO envelopes have so 
far mostly assumed a slab or s pherical s ymmetry and 
a static physical structure (e.g. IViti &: W illiams 199^ 
iDotv et al1l2002n . In situations with chemical timescales 
longer than dynamical timescales, the chemis try needs to 
be co upled with the dynamical evolution (e.g. iDotv et all 
|2006[ ). A particular ex ample for this situation is the 
slow sulfur chemistry (jWakelam et al.l [2004). Other 
models focus on the effect of protostellar FUV and X- 
ray radiation on the chem istry of a spherical envelope 
(jStauber et all 12004 I2005D . Due to the high extinction 
of FUV radiation by dust, they find the chemistry to 
be modified only in the innermost few hundred AU. X- 
rays on the other hand, have a much lower absorption 
cross section. They may thus penetrate further into the 



envelope and dominate the ionization in the entire hot 
core. To study the effects of an outflow cavity that allows 
FUV radiation to escape and irradiate a larger part of 
the envelope, two-dimensional axisymmetric models are 
necessary. 

Non-spheric al chemical models of YSO envelo p es are 
presented by iBrinch et all (|2008D . iVisser et al.] ()2009f ) 
and in paper II. These models assume axisymmetry. The 
first two are applied on low-mass sources and include the 
effects of infall to study the chemical evolution from the 
cloud to the disk. The model in paper II is applied on a 
high-mass star forming region. It calculates the gas tem- 
perature self-consistently with the chemical abundances 
similar to models of photo dominated regions (PD Rs; e.g. 
lKaufmanet~aI1 119991 iMeiierink fc SpaansI [20051) . but it 
does not account for the dynamical evolution. The model 
presented in paper II considers FUV {6 < E^ < 13.6 eV) 
and X-ray {E^ > 100 eV) induced chemistry. 

In this work, we further develop the model introduced 
in paper II. It is based on the grid of chemical mod- 
els introduced in the fi rst paper of this series of papers 
(jBruderer et al.|[2"009d hereafter paper I). We extend the 
model with two important tools: (i) with a dust radia- 
tive transfer method for a self-consistent calculation of 
the dust temperature and (ii) a multi-zone escape prob- 
ability method to study the molecular excitation and ra- 
diation. We apply the model to AFGL 2591 to answer 
the questions raised above. In particular we predict line 
fluxes of two selected diatomic hydrides (SH and CH+), 
having typical properties of diatomic hydrides with high 
temperatures needed for formation and high critical den- 
sities. These two molecules will be observed in AFGL 
2591 by the Water In Star-forming regions with Herschel 
(WISH'^; van Dishoeck et al. 2010, in prep.) guaranteed 
time key program. 

This paper is organized as follows. We first describe the 
details of the model in section [21 The physical structure 
of different models is discussed in section [Hj The chem- 
ical abundances obtained from these physical structures 
are presented in section[4l and compared to the results of 
one-dimensional models. The molecular excitation and 
line radiation is then studied in section [5l and the scien- 
tific prospectives of hydride observations with Herschel 
are discussed in [Sj Limitations of the model used in this 
work are given in Section [Tj 

2. MODEL 

In this section, we briefly describe the calculation flow 
of the models used here. Details on the new parts of 
the model, the dust radiative transfer and the molecular 
excitation calculation, are given in Appendicies \K\ and iBl 
respectively. 

The modeling process starts with assuming a density 
structure and the spectrum of the central source. A 
dust radiative transfer calculation then solves for the 
dust temperature. The local FUV and X-ray radiation 
is calculated at every point of the model as described in 
Appendix [Al In the next step, the gas temperature is ob- 
tained by equating the heating and cooling rates (paper 
II). Since these rates depend on the chemical abundances 
(e.g. of the atomic coolants C"*" and O) the calculation 
of the gas temperature has to be performed simultane- 

^ www.strw.leidcnuniv.nl/WISH 
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ously with the chemistry. This time-consuming step is 
facilitated by the grid of chemical models introduced in 
paper I. The grid consists of a database of precalculated 
abundances depending on physical parameters (density, 
temperature, FUV flux. X-ray flux, cosmic-ray ioniza- 
tion rate and chemical age; see paper I for a definition 
of these parameters) from which abundances can be in- 
terpolated quickly. The interpolation approach is also 
used to obtain the abundances of all other species (e.g. 
diatomic hydrides). The molecular excitation is calcu- 
lated using the escape probability approach described in 
Appendix [B] and yields the modeled line radiation, thus 
the observable quantities. 

2.1. A grid of density structures/cavity shapes 

A spherical model of AFGL 259 1 has been constructed 
bv lvander Taket all (|1999l 120001 ) based on maps of the 
dust continuum, the spectral energy distribution (SED) 
and molecular line emission. In their fitting, they find 
good agreement on larger scales. However, continuum 
emission at shorter wavelengths (warm dust) requires a 
shallower density profile in the innermost 10" — 20". They 
conclude that an inner region with roughly constant tem- 
perature may be requi red to reproduce all observations. 
|j0rgensen et al.l (|2005[ ) encountered similar problems of 
an excess emission in short wavelengths when fitting the 
continuum emission of IRAS 16293-2422. They con- 
cluded that a la rge spherical c a vity c an explain this ex- 
cess. However, iCrimier et al.l (|2010( ) can explain both 
SED and maps by a higher bolometric luminosity. In 
AFGL 2591, new hi g h res olution observations at 24.5 
/im by Ide Wit et al.l (|2009[ ) show extended emission at 
scales of up to 10" in the direction of th e outflow. This 
emissi on stems from the outflow walls (jBruderer et al.l 
l2009a( ). supporting the scenario of a cavity. 

For this work, we adopt the density str ucture of the 
spherical model bv lvan der Tak et al.l(|1999[ ). This choice 
allows to study the differences of the two-dimensional 
chemical model including an o utflow cavity with the 
spherically syrametric mo dels by iDotv et aD (|2002| ) and 
iStauber et al.l(|2004 120051 ). The cavity shape is assumed 
to follow a power-law with index b, 
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where z and r are coordinates along the outflow axis and 
perpendicular to the outflow, respectively. The full open- 
ing angle a of the outflow is given at a distance zq. As in 
paper II, we use zn = 10000 AU , approximately the scale 
on which iPreibisch et al.ll2003l have probed the opening 
angle of the outflow cavity. Cavity shapes for different 
choices of a and b are given in Figure [TJ The solid line 
indicates the model used in paper II. Figure [T] also shows 
the line of sight toward us cro ssing the outflow wall at a 
shallow angle, as suggested bv lvan der Tak et al.l (|1999l ). 
A further parameter of the model is the density in the 
cavity, which we assume to be reduced by a constant 
factor 7 = Uout/'T-in compared to the spherical density 
model at the same radius. This is motivated by the pres- 
sure equilibirum along the outflow wall (paper II). Note 
that for 7 = 2.5 x 10^^, used in most of the models, the 
results are independent of that assumption as the opacity 
in FUV wavelengths along the outflow is low. 
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Fig. 1. — Separation between envelope and outflow using the 
power-law des cription (Equation [T]l . The line of sight suggested by 
Ivan der Tak et al. (1999 ) is given by an arrow. 

A grid of different models is discussed in the follow- 
ing sections. Starting from the model used in paper II 
(referred as model Standard in the following) we explore 
the parameter range by varying the outflow opening an- 
gle a, the power-law index b, cavity density 7 and the 
bolometric luminosity of the protostar Lboi (Table [1}. 
As in paper II, we assume the central source to emit a 
black b ody spectrum with a tem perature of 3 x lO'' K fol- 
lowing l^ndeE^keiil] (|1999f ) (see Section [T]). Thus, 
the FUV field has photons that can dissociate CO and 
H2. 



TABLE 1 

Parameters of the adopted models. 



Model 


a 

n 


b 


7 = "out /"in 


Lhoi 


Standard 


62 


2 


2.5(-3) 


2(4) 


Disk'' 


62 


2 


2.5(-3) 


2(4) 


Opening angle 










a20 


20 


2 


2.5(-3) 


2(4) 


a40 


40 


2 


2.5(-3) 


2(4) 


a 80 


80 


2 


2.5(-3) 


2(4) 


alOO 


100 


2 


2.5(-3) 


2(4) 


Cavity density 










7O 


62 


2 





2(4) 


-fO.Ol 


62 


2 


l(-2) 


2(4) 


-yO.l 


62 


2 


0.1 


2(4) 


Cavity shape 










bl.5 


62 


1.5 


2.5(-3) 


2(4) 


b2.5 


62 


2.5 


2.5(-3) 


2(4) 


Source luminosity 










L2e3 


62 


2 


2.5(-3) 


2(3) 


Lle4 


62 


2 


2.5(-3) 


1(4) 


L4e4 


62 


2 


2.5(-3) 


4(4) 



Note: a{b) = a X lO'' 

° Parameters like in the model Standard, but with a disk of mass 

0.8 Mq in the innermost 400 AU. 



A compact structure with a systematic veloc- 
ity gradient has been detected in AFGL 2591 by 
Ivan der Tak et ahl (|2006D using interferometric observa- 
tions. They propose a circumstellar disk to be the source 
of this radiation and constrain a mass of « 0.8 Mq within 
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a radius of 400 AU to the protostar. While a detailed 
modeling of the disk is beyond the scope of this study, 
model Disk implem.ents a toy-m.odel of a disk em.bedded 
in the envelope of the model Standard. The toy-model 
consists of a massless (in relation to the protostar) Keple- 
rian disk (e.g. Fisc her et al. 1996 ) with a surface density 
E oc r^°-^^, a scale height r/z = 0.15 and a mass of 0.8 
Mq within a radius of 400 AU to the protostar. The 
density distribution of the disk and the inner envelope is 
given in Figure [SJ 

3. PHYSICAL STRUCTURE 

Different physical conditions such as gas and dust tem- 
perature, density, FUV and X-ray irradiation govern the 
chemical composition and molecular excitation. As a 
first step to quantify these physical properties of the dif- 
ferent envelopes models, we study the amount of material 
with distinct conditions (e.g. a gas temperature above 
100 K). Certain conditions like warm temperatures, a 
high density and a FUV field stronger than 1 ISRF (in- 
terstellar radiation field) can be necessary for the forma- 
tion and excitation of a molecule. Combinations of these 
regions (e.g. high density and strong FUV field) can be 
relevant for the formation and excitation of molecules 
and thus the predicted line strength that can be com- 
pared with observations. While the definition of such 
physical regions is arbitrary, it gives a simple overview 
on differences between models and thus the influence of 
the geometry. 

In Table m we show the absolute mass (in Mq) and 
mass relative to Model Standard (Af/Mgtandard) for the 
following regions: 

• TdustjTgas > 100 K. Several molecules only form 
at high temperature, like for example the diatomic 
hydrides discussed in this work. We report on the 
region with dust and gas temperature above 100 
K, which is assumed to be the water evaporation 
temperature (paper I). 

• FUV. FUV radiation alters the chemical compo- 
sition by photoprocesses (ionization and dissocia- 
tion) and also through heating of the gas. The 
FUV region is defined by the FUV flux larger than 
1 ISRF. We recall, however, that the structure of 
the outflow walls is more complicated and some 
species only exist in different layers, for example 
the C+/C/CO layer structure also seen in PDRs. 

• X-rays. X-rays affect the chemistry mostly by 
secondary elec trons created by fast photoelectrons 
ioniz ing H2 ( Malonev et al.l 119961 : iStauber et al.l 
120051 ). Assuming an X-ray luminosity of 10'^^ 
erg s~^ in the 0. 1 - 10 keV band as suggested 
by IStauber et al.1 ()2005[ ). h eating by X-rays ca n 
be neglected (paper II and IStauber et al.1 l2005l ). 
While the influence of X-rays is very similar to 
that of a high cosmic-ray ionization rate (paper 
I), there ar e differences hi FU V irradiated re- 
gions (e.g. IStauber et al.l 120061 for a discussion 
of the water abundance in FUV and X-ray ir- 
radiated regions). The X-ray region is defined 
by the H2 ionization rate larger than the "stan- 
dard" cosmic-ray ionization rate of Cc.r. = 5.6 x 



10 
10 



-17 
-16 



s ^ by a factor of about five (C > 3 x 
s~^). This anticipates the uncertainty in 
the observational determination of the cosmic-ray 
ionization rate ([v an dcr Tak & van Dishocck,, 2000l 
lDotveraIll200a]rndriolo et al.ll2007D . 

• T > 100 K & No FUV. An important region of 
the envelope is the part which has a temperature 
above 100 K but no FUV irradiation. Hot core 
species which can be photodissociated may exist 
only in that region. The region is defined by the 
temperature T = max(Tdust,T'gas) (see Section [5]) 
and the FUV region, defined above. 

The spherical model is also given in Tabl e |21 For com- 
parison with earlier modeling, we follow IStauber et al.l 
(|2004 120051 ) and assume a visual extinction of ry = 
and Go = 10 ISRF at the inner edge of the modeled re- 
gion at 200 AU an d use the gas and dust temperature by 
iDotv et al.1 (|200l . 

For the model Spherical, a plot of the radial depen- 
dence of density and temperature is given in Figure 3 
of paper I. For this model, the radial size of the regions 
discussed here are: 1900 AU (Tdust > 100 K), 1100 AU 
(Tgas > 100 X and T > 100 K (No FUV)), 224 AU 
(FUV), 2300 AU {X-rays and only X-rays). This shows 
the small size of the FUV irradiated and heated region 
in a spherical model compared to the models including 
and outflow wall. The extent of the different regions of 
model Standard is indicated in Figure [2] and the inset of 
Figure [6l 

The X-ray luminosity of AFGL 2591 cannot be mea- 
sured directly do to the high column density towards the 
pr otostar, but upper lim its of 1.6 x 10"^^ erg s~^ are given 
in lCarkner et al.l(|1998f) . However, they assume an atten- 
uating column density to the source of only iVn = 10^^ 
cm~^, much smaller than the radial column density of 
about iVn = 3 x lO^'' cm^^ of the density proflle used 
here. Note that the column density derived from the den- 
sity proflle can be aff ected by uncertainties e .g. in the 
adopted dust opacity (jvan der Tak et al"1ll999n . For the 
plasma temperature of 7x lO'' K assumed here (paper II), 
the X-ray luminosity may be a factor of about 10 higher 
assuming an attenuating column density of TVh = 3 x 10^'^ 
cm^^ (e.g. Appendix A of paper I). We thus consider the 
adopted value of the X-ray luminosity of 10^^ erg s~^ as 
an upper limit. 

3.1. Total mass 

The total mass of the envelopes (Aftot) of the various 
models is very similar with maximum differences of 25 
% compared to model Standard (Table [2|) . The situa- 
tion is different if only the inner part of the envelope is 
considered. We give the mass within a radius of 4000 
and 20000 AU to the protostar (M4000 and M20000) fol- 
lowing the Herschel beam at the shortest/longest wave- 
lengths of PACS/HIFI (60-625 p,m.) assuming a distance 
of 1 kpc to the source (|van der Tak et 311 11999'). Clearly, 
there is foreground material within the Herschel beam. 
In this work, however, we concentrate on molecules with 
lines having high critical densities. They thus trace only 
the inner regions due to the density gradient. However, 
dust or line radiation from inner regi ons may also pu mp 
molecules further out (e.g. ,van Kempen et al.ll2008f) . A 
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TABLE 2 

Mass of different regions in the 2D models and the spherical model. The upper part of the table gives absolute values of 

THE mass (in Mq), while THE LOWER PART OF THE TABLE GIVES FACTORIAL DEVIATIONS " RELATIVE TO THE MODEL Standard. FACTORIAL 

DEVIATIONS LARGER THAN A FACTOR OF THREE ARE HIGHLIGHTED. THE TOTAL MASS OF THE ENVELOPE AND THE MASS WITHIN A RADIUS OF 

4000 AU AND 20000 AU to the PROTOSTAR are given by Mtot, A^4000 and M20OOO. respectively. "1" MEANS THAT THE VALUE IS BY 

DEFINITION EQUAL TO THE MODEL Standard. THE SPATIAL EXTENT OF THE REGIONS IS DISCUSSED IN SECTION [3] 



Model Mtot 


A'^4000 


A^20000 


Tduet > 100 K 


Tgas > 100 K 


FUV 


X-rays 


only X-rays'' 


T > 100 K (No FUV)'= 


Absolute Values (Mq) 
















Spherical 49.4 


0.94 


33.5 


0.20 


6.8(-2) 


5.7(-4) 


0.29 


0.29 


6.7(-2) 


Standard 44.3 


0.45 


18.7 


6.7(-2) 


0.71 


1.9 


0.70 


0.32 


7.9(-3) 


Relative to Model 


Standard 
















Spherical 1 . 1 


2.1 


1.8 


3.0 


9.6(-2) 


3.0(-4) 


0.42 


0.92 


8.5 


Standard 1 


1 


1 


1 


1 


1 


1 


1 


1 


Disk 1.1 


9.1 


1.2 


15.6 


2.3 


1.2 


1.3 


0.82 


83.7 


Opening angle 


















a20 1.1 


1.8 


1.2 


2.3 


0.36 


0.25 


0.70 


1.3 


17.0 


a40 1.1 


1.5 


1.1 


1.5 


0.60 


0.60 


0.87 


1.2 


6.7 


a80 0.91 


0.64 


0.87 


0.72 


1.4 


1.4 


0.90 


0.47 





alOO 0.75 


0.36 


0.68 


0.49 


1.8 


1.8 


0.73 


9.4(-2) 





Cavity density 


















■yO 1 


1 


1 


1.1 


1.2 


1.1 


1.0 


1.0 


0.71 


•yO.Ol 1 


1 


1 


0.91 


0.46 


0.73 


0.93 


1.0 


1.8 


■yO.l 1 


1 


1 


0.71 


1.2(-2) 


3.8(-3) 


0.65 


1.4 


5.5 


Cavity shape 


















bl.5 0.99 


1.3 


1.01 


1.2 


0.90 


0.95 


0.76 


0.85 


4.3 


b2.5 1 


0.76 


0.99 


0.91 


1.1 


1.1 


1.1 


0.90 


8.9(-2) 


Source luminosity 


















L2e3 1 


1 


1 


6.9(-2) 


0.37 


0.66 


1 


1.2 





Lle4 1 


1 


1 


0.41 


0.76 


0.89 


1 


1.0 


8.9(-2) 


L4e4 1 


1 


1 


2.3 


1.3 


1.1 


1 


0.96 


5.0 



Note: a{b) = a X 10*. 

"Defined by niaa;(M/Mstandard, A'-fstandard/^) > 3. 

''Only in the X-ray region and not in the FUV region. 

'^Region with T = max(Tgas, Tdust) > 100 K outside the FUV region. 



detailed discussion which regions are traced by molecu- 
lar emission requires the calculation of the excitation as 
presented in Section [5j 

Models with a large opening angle (a 100) or the toy- 
disk have masses M4000 which are 3 times smaller and 10 
times larger, respectively, compared to the model Stan- 
dard. More than half of the envelope mass is outside a 
radius of 20000 AU and differences between individual 
models are less evident in M20000 and not seen in Mtot- 

We conclude that different cavity shapes affect the to- 
tal mass only in the inner part of the envelope, with the 
largest changes found for models with a large/small out- 
flow opening angle a or a disk embedded in the envelope. 

3.2. Dust temperature above 100 K 

The mass with dust temperature above 100 K approx- 
imately scales with M4000 for models with the disk and 
different cavity shape (Models a20-al00, hi. 5, b2.5). 
This is explained by the dust temperature depending 
only on the distance to the protostar and not on the 
density in optically thin, spherically symmetric models. 
For this c ase, the radius for a giv en temperature varies as 
-^boi (^-S- lJ0rgensen et al.|[2006[ ). For the density struc- 
ture oc r~^, the mass with Tdust > 100 K should scale 
oc Lboi- Models L2e3, Lle4 and L4e4 depend slightly 
stronger on Lboi due to the cavity removing more of the 
dense material in the inner part compared to larger dis- 
tances to the protostar. Changing the density in the cav- 
ity (Models 7O, J 0.01 and ■'yO.l) affects the mass with 



Tdust > 100 K by the attenuation of direct stellar ra- 
diation in FUV wavelengths that heats the dust in the 
outflow walls. 

We conclude that the amount of material with dust 
temperature above 100 K is within a factor of two to the 
model Standard for models with a different cavity shape, 
but larger for the models including a Disk and with a 
lower bolometric luminosity (L2e3). 

3.3. Gas temperature above 100 K 

Gas with Tgas > 100 K is more than an order of mag- 
nitude less present in the spherical model and model 
^0.1 compared to the model Standard. This underlines 
the impact of the empty cavity to extensively enhance 
warm gas by escaping FUV radiation. Model ^0.1 has 
even less warm gas than the spherical model, due to a 
more efhcient cooling through the low-density outflow re- 
gion. In all other models, the deviation of the mass with 
T'gas > 100 K compared to model Standard is smaller 
than a factor of three. In particular, the models with a 
different bolometric luminosity {L2e3, Lle4 and L4e4) 
scale less than with Lboi due to the PDR surface temper- 
ature scaling wit h less than Go fo r densities be low about 
lO^cm"^ (e.g. Kaufman et all [T999 . Mciierink et al.l 
12001 . 

We conclude that the mass with gas temperatures 
above 100 K changes most if the geometrical situation 
is similar to the spherically symmetric model. These are 
models with a small outflow opening angle or a high den- 
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sity of absorbing material in the outflow. 

3.4. X-ray and FUV irradiation 

Regions with high-energy irradiation (FUV and/or X- 
rays) are given in Figure [5] for models Standard, disk and 
a 20. They show an FUV region along the outflow cav- 
ity. The X-ray region is less concentrated to the surface 
and also includes the envelope to distances of up to a few 
1000 AU. A region with only X-ray influence follows the 
FUV irradiated surface layer in the vicinity of the proto- 
star (distance < 5000 AU). Some surface concentration 
is also seen in X-rays due to photons with energy be- 
tween 0.1 and 1 keV which are attenuated by a relatively 
low column density (of order 10^^ cm~^). Harder X-rays 
with Ej > 1 keV require a column density of order 10^"* 
cm~^ to be attenuated (paper I). A disk can provide such 
a high column density and shield the envelope even for 
hard X-rays (Figure [21). 



a.) Standard (inner region) 



b.) Disk 



f 



2 ■ 10"' 



I 



1-10* 




Fig. 2. — FUV and X-ray irradiated regions in the model Stan- 
dard and the models with a disk and a = 20°. FUV and X-ray 
irradiated regions as defined in Section [3] are given in green (only 
X-rays), orange (only FUV) and red (both X-rays and FUV). 

The mass of the X-ray/FUV irradiated gas is the same 
for most models within a factor of two relative to model 
Standard. Exceptions are the FUV irradiated mass in the 
spherical model and model jO.l being orders of magni- 
tude smaller. Smaller opening angles (Model a20) yield 
a larger incident angle of the FUV radiation and thus a 
thinner outflow wall. This is reflected in the FUV irradi- 
ated mass and also the mass with Tgas > 100 K. For large 
opening angles {a 100) the FUV radiation penetrates to 
the midplane of the flat structure and the region with 
only X-ray irradiation is much smaller. In the model in- 
cluding a disk, the mass suffering X-ray irradiation is not 
much different despite the shielding seen in Figure [2] It 
is caused by the larger density in the innermost part at 
the surface of the disk. 

The largest variations between different models are 
found for the mass of the region with T > 100 K but 



no FUV irradiation. This region is shown in yellow in 
the Figures [6] and [5] for models Standard, hi. 5, h2.5, a20 
and disk. A thin layer of FUV heated gas is followed by 
an infrared heated region. If the region close to the pro- 
tostar is geometrically thin {b2.5, a80 and alOO) FUV 
radiation can penetrate to the midplane and no material 
with T > 100 K but no FUV irradiation exists. Geo- 
metrically thin/thick means that the distance in verti- 
cal direction between the midplane and the outflow is 
small/large. In models with a geometrically thick inner 
region {hi. 5, a 80 and a 100) or more mass in the in- 
nermost region [disk) the amount of IR heated gas with 
T > 100 can be even larger than in the spherically sym- 
metric model. The strong sensitivity of this mass on the 
bolometric luminosity is the result of the dust mass with 
temperature above 100 K depending more on the bolo- 
metric luminosity than the gas temperature (Section 13.21 
and ESI). 

We conclude that the mass with X-ray irradiation is 
not considerably affected by the geometry, while the mass 
with FUV irradiation shows differences for those models 
that are similar to the spherical situation (small open- 
ing angle or high density of absorbing material in the 
outflow). The region with temperature above 100 K but 
no FUV irradiation depends considerably on the geome- 
try and the bolometric luminosity of the source. Models 
with a geometrically thin inner part of the envelope (for 
example a large outflow opening angle) have consider- 
ably less material with temperatures above 100 K but no 
FUV irradiation. 

4. CHEMICAL ABUNDANCES 

Chemical abundances, their evolution, differences be- 
tween physical models and the impact of X-rays are 
studied in this section. We first study fractional abun- 
dances n{y^,r')/ntot{r'), with n{y^,r') and ntotir'), the 
density of the species and the total hydrogen density 
[ntot = 2ri(H2) -l- n(K)) depending on the position. The 
fractional abundance gives an overview of the spatial 
variation of the molecule. Next, the volume averaged 
fractional abundance is discussed. The volume averaged 
abundance of a species X within a radius r is defined by 



{X)r 



n(X,r')dV" 



\r' I <r 



' \r' I <r 



ntot{r')dV . (2) 



This quantity is more related to observations, as the vol- 
ume averaged abundance times the mass gives the total 
amount of molecules. The total amount of molecules 
scales with the line flux in the optically thin case, ne- 
glecting excitation effects. 

In paper I, the dust and gas temperature are assumed 
to be equal, and evaporation is approximated by an in- 
creased initial abundance for certain temperature ranges. 
Here, we calculate the dust and gas temperature ex- 
plicitly. Since the abundances of the species consid- 
ered here do not change much in the temperature range 
70 < T < 100 K, but much more with water evapora- 
tion at 100 K, we use T — max(T'dust,'?gas) to obtain 
the abundances from the grid of chemical models. This 
choice affects only the innermost region with dust tem- 
peratures above 100 K heated by infrared photons (Fig- 
ure |3] and Section [3]). Abundances outside this region 
are found to agree to within 20% if the gas temperature 
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were used instead. This also shows that the abundances 
of the diatomic hydrides are not sensitively depending 
on the assumed water evaporation temperature. The 
evaporation temperature of the diatomic hydrides is not 
well known. Howeve r, as their ab undance on du st grains 
likely is small (e.g. iHasegawa fc Herbstlll993l ). the ef- 
fect of direct evaporation of diatomic hydrides to the gas 
phase is neglected here. 

4.1. Standard and spherically symmetric model 

The abundances of diatomic hydrides, related species 
and main molecules/atoms (H2, CO, CO2, H2O, O, C, 
C"*") for the model Standard are given along cuts of con- 
stant z (parallel to the midplane, at 0, 2000 and 10000 
AU) in Figures |3] and |4l Presented abundances are for 
a chemical age of 5 x 10^ yrs (jStauber et al.l [20051) and 
with a protostellar X-ray luminosity of 10'^^ erg s~^ or no 
protostellar X-rays. The density structure of the model 
along with the dust and gas temperature, total ioniza- 
tion rate (see paper I), the FUV flux and attenuation 
are given in the top panels of Figure |3l 

As a rough classification of the abundance profiles we 
can divide the species into broad groups of molecules 
which are (i) enhanced in the outflow wall compared to 
the envelopes (e.g. electrons, C"*", CH"*") with different 
widths of the enhanced region; (ii) destroyed along the 
outflow (e.g. H2, CO or H2O); (iii) enhanced in a thin top 
layer and destroyed at intermediate depths (e.g. OH or 
HaO"*"). Why are some molecules enhanced or destroyed 
in different layers? The overall chemistry of FUV/X-ray 
irradi ated molecular gas has been studied in various pa- 
pers ([HoUenbach fc TielensI 119991 f or a review) and for 
envel opes of YSOs in particular by iStauber et al.l ()2004 
120051 ) . Here, we focus on the chemistry of diatomic hy- 
drides and related species in the physical regime of the 
outflow walls. 

4.1.1. The chemistry in the outflow walls 

In the top layer of the outflow wall with Ty smaller 
than a few, the temperature is high, with Tgas > 100 K. 
Due to self-shielding, H2 and CO are only destroyed in 
a thin layer. Photoprocesses proceeding through contin- 
uum absorption, as for example the photoionization of 
carbon, require a minimum photon energy (1102A). Due 
to the higher dust opacity at shorter wavelengths (red- 
dening), photons with higher energy are absorbed more 
quickly and carbon photoionization is most active at the 
surface. Ionized carbon forms CH+ by, 0+ -I- H2 — >■ CH+ 
-f- H, but this reaction is highly endoergic by 4640 K. 
The reaction of CH+ with H2 produces CH^ and CH;^ 
which recombine with electrons to CH. Since the oxygen 
ionization energy is slightly higher than 13.6 eV, 0+ is 
produced by charge exchange with protons, O + H+ — >■ 
0+ + H, with H+ from the photodissociation of CH+. 
Ionized hydroxyl (0H+) is then formed by H2 -1-0+ — ?> 
0H+ -I- H. Subsequent reactions with H2 lead from 0H+ 
to H2O+ and then H3O+, which recombines to OH and 
H2O. Hydroxyl is however mainly formed by the endoer- 
gic reaction H2 + O + 3150 K ->■ OH -I- H. Water also 
forms by an endoergic reaction, H2 -I- OH -I- 1740 K — ;> 
H2O -I- H and is photodissociated to OH + H. Reactions 
of OH with O and CO form O2 and CO2, respectively, 
which are both quickly photodissociated. Sulfur hydride 



and its ion are formed in the endoergic reactions S -I- H2 
+ 6130 K ^ SH + H and S+ -F H2 -H 9860 K -^ SH+ 
-\- H with S+ from direct photoionization of S. Finally, 
NH is created by H2 -I- N -I- 18100 K ^ NH + H. Since 
the ionization potential of atomic nitrogen is higher than 
13.6 eV, NH+ is only produced by photoionization of NH 
or charge exchange H+ -1- NH — > NH+ + H. 

In an intermediate layer, at higher optical depth, the 
temperature drops below 100 K (Figure[3]). In this colder 
environment, the photodissociation of OH is still consid- 
erable but the endoergic formation reaction of H2 -1-0 
-> OH -I- H is no longer efficient. This explains the drop 
of the OH abundance in this layer. The drop is also seen 
in other molecules (SH+, NH, NH+, H2O; group (iii) in 
Section |4TT|) . For H3O+, the situation is different: Since 
CH+ is not available in this layer, the main production 
of H+ is disabled and subsequently, no 0+ is available. 
Thus, the abundances of 0H+, H2O+ and H3O+ are de- 
creased. 

We conclude that the concentration of many outfiow 
wall enhanced species towards the outflow is a conse- 
quence of (atomic) ionization being only active at the 
surface due to a quick shielding of short wavelength pho- 
tons by dust (0+ and S+) combined with the high re- 
quired temperature to overcome the activation energy 
for the formation (especially CH+, OH, SH, SH+, NH 
and NH+). 

In the midplane (z = 0), at a distance of about 1000 
AU to the protostar, the influence of H2S, H2O and CO2 
evaporation can be seen. At this point, FUV radiation is 
present but too weak to heat the gas above 100 K, while 
IR radiation still heats the dust grains to Tdust > 100 K. 
The evaporated water is quickly photodissociated and its 
abundance is not increased. The higher amount of sul- 
phur in the gas phase leads to an enhancement of SH and 
SH+ which can been seen. Also affected are C+ and CH, 
which have an increased abundance in colder regions due 
to reaction partners being less abundant. 0+ is destroyed 
in reactions with atomic sulphur, while CH is destroyed 
with atomic oxygen. Among the species discussed here 
(SH, SH+, CH and CH+), only SH reaches abundances 
in this region that are higher compared to the outflow 
walls. Thus, SH emission may be dominated by this re- 
gion wi th su lphur evaporated to the gas phase (Section 
14:2] and lO)) . 

We next consider HCO"*" which is p redicted to be an 
io nization tracer fe.g. | van der Tak fc v an Dishoeck 200Q 
or iSavage fc Ziurvsl l2004| ) . In the top layer of the out- 
flow wall, it is formed by the reaction of II2 with HOC+ 
and C0+. The huge enhancement of C0+ by four or- 
ders of magnitude in a thin layer along the outflow has 
been discussed in paper II. Going deeper into the outflow- 
wall, HCO+ is destroyed by electron recombination and 
then enhanced by the reaction of CO with H;^. The 
Hj" ion is a key spec ies for the X-ray drive n chemistry 
(jMalonev et al.lll996l or lStauber et al.ll2005D and formed 
by cosmic-ray or X-ray ionization of H2 followed by a 
quick reaction with H2. The high abundance of H3 at 
the surface of the outflow wall is maintained by the re- 



action NH+ + H2 



N 



H3 . Deeper in the outflow 



wall, the high electron fraction leads to quick electron 
recombination of H^ . 
An important difference between X-ray and FUV 
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Fig. 3. — Abundances along cuts of constant z through the model Standard. Abundances assuming no protostellar X-ray radiation (solid 
line) and an X-ray luminosity of 10''^ erg s~^ (dashed line) are shown (bottom panels). The density structure and dust/gas temperature, 
total ionization rate and FUV flux/attenuation along cuts with z = 0, 2(300, 10000 AU are given too (top panels). The standard cosmic-ray 
ionization rate of 5.6 X 10~ s~ is indicated by a dashed line. 



driven chemistry is the electron fraction. For fluxes of 
1 erg s^i cm-2 (Gq « 1000 ISRF) the electron frac- 
tion at a total density of 10^ cm^'^ is approximately 
x{e-) = n{e-)/ntot '-' 10"'' for FUV but only IQ-^ for 
X-rays. Since the luminosity in the FUV band is 4 orders 
of magnitude higher than the X-ray luminosity, the large 
electron fraction is mainly produced by FUV radiation. 
Since this electron fraction leads to the destruction of 
H;^ , the influence of X-rays is dwarfed by FUV radiation. 
In ad dition, many molecule s predicted to be X-ray trac- 
ers bv lStauber et ahl ()2005| ). are much more enhanced by 



FUV compared to X-rays. Particular examples are CH+ 
and OH"*", which both have fractional abundances below 
10~^^ in the spherical model but increase to more than 
10~* in the outflow walls. We conclude that the effect of 
X-rays compared to FUV on the chemistry is small in a 
geometry which allows the much stronger FUV radiation 
to escape from the innermost part. 

4.1.2. Volume averaged abundances 

How is the total amount of a molecule/atom altered 
by the enhancement or absence in the outflow wall? For 
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Fig. 4. — Abundances along cuts of constant z through the model Standard. Abundances assuming no protostellar X-ray radiation (solid 
line) and an X-ray luminosity of 10"^^ erg s^^ (dashed line) are shown. 



example a thin layer of strongly enhanced C0+ is found 
in paper II to enhance the total amount by three or- 
ders of magnitude in agreement with observations. We 
note that C0+ is related to the diatomic hydrides stud- 
ied here through formation by C+ + OH — j- C0+ + H. 
FUV destroyed species on the other hand, should not 
decrease their total abundance significantly since only 
a relatively small fraction of the envelope mass is cut 
out by a cavity (Section [3]) . In Table [3l we give vol- 
ume averaged abundances within a radius of 4000 and 
20000 AU for the model Standard and compare it to the 
spherical model. In the following we consider a deviation 
larger than a factor of three to be significant, since simi- 
lar models (jDotv et al.ll2004| ) qualify a factorial deviation 
of this factor between a model and observation as a good 
agreement. 

Outflow wall enhanced species indeed have a signifi- 
cantly enhanced volume averaged abundance compared 
to the spherically symmetric model. Most enhanced are 
the ionized molecules C+, CH+, 0H+, NH+, H2O+ and 
the electron fraction with an enhancement larger than 3 
orders of magnitude and up to almost 7 orders of mag- 
nitude. On the other hand, it is found that the FUV 



destruction of H2O, CO2 and SH reduces the volume av- 
eraged abundance less, by at most a factor of 100. Both 
enhancement and reduction are seen more in the abun- 
dance averaged over 4000 AU rather than 20000 AU. The 
width of the region with FUV enhancement is important 
too. For example the volume averaged abundance of SII+ 
is only increased by a factor of about two despite the en- 
hancement by about 2 orders of magnitude in a very thin 
layer on the surface of the outfiow wall. Similarly, the 
abundance of NH is only increased in the inner part of 
the envelope due to the strong temperature dependence 
of the abundance and the volume averaged abundance is 
only increased by a factor of two. This strong tempera- 
ture dependence and consequently the concentration to 
the inner part also explains the largest difference of en- 
hancement of NII+ between the volume averages of 4000 
AU and 20000 AU. 

The temporal evolution is studied in Tablc[3]by the vol- 
ume averaged abundances within 20000 AU for a chem- 
ical age of 10" s (3 x 10^ yrs) and lO^^ s (3 x 10^ yrs), 
and comp ared to the chemical age of 5 x 10^ yrs found 
bv Stauber et al.l (p005.) . The abundances averaged over 
4000 AU are not shown since they are less time depen- 
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TABLE 3 

Volume averaged abundances (X),- of selected species for the model Standard ("2D") and the spherically symmetric model 
("ID"). Abundances averaged over a radius of 4000 and 20000 AU from the protostar and assuming different chemical 

AGES AND PROTOSTELLAR X-RAY LUMINOSITIES ARE GIVEN. ThE RATIO OF AVERAGED ABUNDANCES BETWEEN THE MODEL Standard AND 

THE SPHERICALLY SYMMETRIC MODEL IS GIVEN ON THE RIGHT. 



Model 


2D 


2D 


2D 


2D 


2D 


ID 


ID 


ID 


Ratio 2D/ 


ID 


X-rays [erg s^-"-] 


1(32) 


1(32) 


- 


1(32) 


1(32) 


1(32) 


1(32) 


a 


1(32) 


1(32) 


- 


Chem. age [yrs] 


5(4) 


5(4) 


5(4) 


3(3) 


3(5) 


5(4) 


5(4) 


5(4) 


5(4) 


5(4) 


5(4) 


Averaged [AU] 


20000 


4000 


4000 


20000 


20000 


20000 


4000 


4000 


20000 


4000 


4000 


H2 


5(-l) 


5(-l) 


5(-l) 


5(-l) 


5(-l) 


5(-l) 


5(-l) 


5(-l) 


1.0 


0.9 


0.9 


e~ 


8(-6) 


5(-5) 


5(-5) 


6(-6) 


8(-6) 


l(-8) 


2(-8) 


6(-9) 


8(2) 


3(3) 


8(3) 


H^ 


3(-10) 


3(-10) 


5(-ll) 


3(-10) 


3(-10) 


3(-10) 


4(-10) 


6(-ll) 


1.0 


0.7 


0.8 


CH 


6(-10) 


2(-9) 


2(-9) 


9(-10) 


6(-10) 


5(-ll) 


8(-ll) 


2(-ll) 


11 


30 


1(2) 


CH+ 


4(-10) 


4(-9) 


4(-9) 


4(-10) 


4(-10) 


4(-16) 


6(-16) 


6(-17) 


1(6) 


8(6) 


7(7) 


OH 


2(-8) 


l(-7) 


l(-7) 


2(-8) 


2(-8) 


9(-9) 


l(-8) 


3(-9) 


1.9 


9.5 


43 


OH+ 


I(-IO) 


3(-9) 


3(-9) 


I(-IO) 


I(-IO) 


2(-14) 


3(-14) 


4(-15) 


6(3) 


9(4) 


7(5) 


SH 


3(-ll) 


2(-ll) 


2(-ll) 


2(-10) 


2(-ll) 


2(-10) 


l(-9) 


3(-9) 


0.1 


0.01 


6(-3) 


SH+ 


l(-ll) 


2(-10) 


2(-10) 


l(-ll) 


l(-ll) 


9(-12) 


9(-ll) 


2(-12) 


1.1 


2.1 


1(2) 


NH 


l(-9) 


2(-9) 


9(-10) 


l(-9) 


l(-9) 


5(-10) 


7(-10) 


6(-ll) 


2.4 


2.3 


15 


NH+ 


1(-12) 


3(-ll) 


3(-ll) 


1(-12) 


1(-12) 


1(-15) 


2(-15) 


3(-16) 


8(2) 


2(4) 


1(5) 


C+ 


7(-6) 


3(-5) 


3(-5) 


4(-6) 


7(-6) 


4(-10) 


7(-10) 


I(-IO) 


2(4) 


4(4) 


2(5) 


C 


5(-7) 


2(-6) 


7(-7) 


l(-7) 


6(-7) 


4(-8) 


3(-8) 


l(-8) 


11 


70 


51 


CO 


2(-4) 


l(-4) 


l(-4) 


2(-4) 


2(-4) 


2(-4) 


2(-4) 


2(-4) 


1.0 


0.8 


0.8 


C02 


5(-9) 


4(-9) 


5(-9) 


l(-9) 


l(-8) 


3(-8) 


7(-8) 


9(-7) 


0.2 


0.06 


5(-3) 


H20 


2(-7) 


9(-8) 


6(-8) 


5(-8) 


3(-7) 


8(-7) 


6(-7) 


5(-6) 


0.3 


0.1 


0.01 


H20+ 


5(-ll) 


5(-10) 


5(-10) 


5(-ll) 


5(-ll) 


4(-14) 


7(-14) 


1(-14) 


1(3) 


7(3) 


5(4) 


H30+ 


4(-10) 


l(-9) 


9(-10) 


3(-10) 


4(-10) 


l(-9) 


2(-9) 


6(-10) 


0.4 


0.6 


1.5 


HCO+ 


2(-9) 


l(-9) 


4(-10) 


2(-9) 


2(-9) 


4(-9) 


8(-9) 


2(-9) 


0.5 


0.2 


0.2 





5(-5) 


8(-5) 


8(-5) 


5(-5) 


4(-5) 


4(-5) 


3(-5) 


3(-5) 


1.3 


2.5 


2.4 



Note: a(6) = a X lO*". 
"FUV also switched off. 



dent. The outflow wall enhaneed species have very short 
chemical time-scale for their formation and thus do not 
show any temporal evolution. We find that CH+ in the 
enhanced region only requires of order 10 yrs to reach the 
final fractional abundance to within a factor of two. The 
time-scales of 0H+, NH+ and H20''' are similarly short. 
Assuming an upper limit of 10 km s~^ for the Alfven 
and sound velocity, the irradiated gas would move up to 
a distance of about 25 AU within the chemical timescale 
which justifies our approach to ignore the dynamical evo- 
lution of the gas. This is further justified by the fact that 
diatomic hydrides produced in the outflow walls diffus- 
ing into the cloud are quickly destroyed by reactions with 
other species. To test this, a chemical model presented 
in paper I has been run for conditions of the cold and not 
irradiated gas but with initially increased abundances of 
the diatomic hydrides. Since the C^ enhanced region 
is thicker, contributions from colder and less irradiated 
regions with longer timescales are important and the vol- 
ume averaged abundance shows some temporal evolution 
(within a factor of two) . 

4.2. Comparing different geometries 

The influence of different physical structures (Section 
[3]) on chemical abundances is studied in this section. Ta- 
bles m and [5] give the volume averaged abundances of 
molecules studied in the previous section. A chemical 
age of 5 X 10* yrs and an X-ray luminosity of Lx — 10'^^ 
erg s~^ are used. The abundance is averaged over 4000 
AU (aU models) and 20000 AU (selected models). To 
obtain the total amount of a molecule, the mass of the 
region over which is averaged (M4000 and M20000) is 



given relative to the model Standard. The molecules 
which were found in Section 14.1.21 to be signiflcantly en- 
hanced/destroyed or unchanged are marked with "+"/"- 
" or "0", respectively. 

For the outflow wall enhanced species CH+, 0H+, 
SH+, NH"'', C"*", H2O"'" and electrons, the volume aver- 
aged abundance is most changed in the models Disk, a20 
and 7O. i. While for a20 and ^0.1 averaging over 4000 
and 20000 AU leads to a signiflcantly lower abundance, 
the differences in the Disk model are only seen for the 
average over 4000 AU. However, when scaling with the 
higher mass of the Disk model within 4000 AU, the total 
amount of these charged species remains approximately 
constant. This reflects that the fractional abundance of 
outflow wall enhanced species is approximately inversely 
proportional to the density (paper I). In addition, the 
higher density yields a lower gas temperature and thus a 
decreased amount of molecules which require a high tem- 
perature for formation (e.g. CH"*"). Quantitatively, we 
flnd that the volume averaged abundances of outflow wall 
enhanced species behave similarly as the FUV irradiated 
mass, discussed in Section [31 

The FUV destroyed molecules SH, H2O and CO2 have 
a similar dependence of the volume averaged abundance. 
Like the region without FUV irradiation but tempera- 
ture above 100 K (Section [3]), they are more dependent 
on the opening angle {a20-al00) and the shape {hi. 5 
and h2.5) than the FUV enhanced molecules. Quantita- 
tively, however, the dependence of the volume averaged 
abundance on the physical models differs more between 
these three species than for the outflow wall enhanced 
species. This is due to different effects entering the abun- 
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TABLE 4 
Ratios of the volume averaged abundances (Eq. [2J to the model Standard. A chemical age of 5 x 10* yrs and a 

PROTOSTELLAR X-RAY LUMINOSITY OF lO"^^ ERG S"-"" ARE ASSUMED. THE MASS OVER WHICH IS AVERAGED (M400O AND Af2000o) IS GIVEN 

RELATIVE TO MODEL Standard. 



Model 


H2 


e 


H,t 


CH 


CH+ 


OH 


OH+ 


SH 


SH+ 


NH 


Within 


4000 AU 


A^4000 


0" 


+ 





(+) 


+ 


(+) 


+ 


- 


+ 





Spher 


cal 


2.1 


1.1 


4(-4) 


1.4 


3(-2) 


l(-7) 


0.1 


l(-5) 


1(2) 


0.5 


0.4 


Standard 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Disk 




9.1 


1.1 


1(-1) 


3(-2) 


0.2 


7(-2) 


0.7 


6(-2) 


2.5 


8(-2) 


5(-2) 


a20 




1.8 


1.0 


0.2 


1.1 


0.2 


0.2 


0.9 


0.3 


1(2) 


0.3 


1.1 


a40 




1.5 


1.0 


0.5 


1.2 


0.4 


0.5 


0.5 


0.6 


4.0 


0.6 


1.1 


a80 




0.64 


1.0 


1.5 


0.7 


2.9 


1.8 


1.9 


1.2 


0.4 


1.6 


1.0 


alOO 




0.36 


0.9 


2.2 


0.4 


8.0 


2.9 


3.8 


1.5 


0.1 


2.3 


1.2 


7O 






1.0 


1.0 


1.0 


0.9 


1.0 


0.8 


0.9 


0.9 


0.8 


0.9 


■yO.Ol 






1.0 


0.8 


1.0 


0.9 


0.8 


0.8 


0.5 


1.4 


0.9 


1.0 


-yO.l 






1.1 


4(-2) 


1.6 


9(-2) 


9(-3) 


0.2 


l(-2) 


4.0 


0.3 


1.8 


bl.5 




1.3 


1.0 


0.6 


1.0 


0.6 


0.6 


0.8 


0.6 


3(1) 


0.7 


1.4 


b2.5 




0.76 


1.0 


1.4 


1.0 


1.7 


1.6 


1.4 


1.4 


0.6 


1.4 


1.0 


L2e3 






1.0 


0.5 


1.2 


1.7 


0.5 


1.1 


0.1 


0.8 


0.6 


1.2 


Lle4 






1.0 


0.8 


1.1 


1.3 


0.9 


1.1 


0.8 


0.7 


1.0 


1.3 


L4e4 






1.0 


1.3 


0.9 


0.9 


1.2 


0.8 


1.0 


2.7 


0.9 


0.7 


Within 


20000 AU 


Aif20000 






















Spherical 


1.8 


1.0 


l(-3) 


1.0 


9(-2) 


9(-7) 


0.5 


2(-4) 


6.8 


0.9 


0.4 


Standard 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Disk 




1.2 


1.0 


0.8 


0.8 


0.9 


0.7 


1.6 


0.5 


1.0 


0.6 


0.8 


a20 




1.2 


1.0 


0.2 


1.0 


0.4 


0.2 


0.8 


0.3 


5.7 


0.6 


0.7 


-yO.l 






1.0 


l(-2) 


1.1 


0.3 


2(-3) 


0.7 


8(-3) 


1.2 


0.4 


1.1 


L2e3 






1.0 


0.5 


1.0 


1.2 


0.2 


0.8 


8(-2) 


1.3 


0.5 


1.0 



Note: a{b) = a X 10*. 

"Ratio to the spherical model within 4000 AU (Table [S] + more than lO'', 



less than 0.1, (+) more than 3, within 0.1 to 3). 



dances of SH, H2O and CO2, like evaporation combined 
with a different depth dependence of the photodissocia- 
tion rates. For example, the larger dependence of CO2 
on the opening angle for models a80 and a 100 is a result 
of the concentration of the molecule towards the center 
combined with the CO2 photodissociation rate dropping 
faster with extinction compared to H2O and SH. On the 
other hand, SH is more enhanced in the models a20 and 
bl.5 due to the factor of 10 larger increase in abundance 
in the hot-core region. Quantitatively, photodissociated 
and hot-core enhanced species thus have a more com- 
plex dependence on the geometry than FUV enhanced 
species. 

Shocks, which are not accounted for here, may release 
H2O and CO2 to the gas phase. In a fast J-shock with 
vs = 80 km s~^ and a density of 10^ cm~'^, the width of 
dust grain s with temperature ab ove 100 K is of order 700 
AU (.HoUenbach fc McKedll989f ). Thus along the outflow 
walls, H2O and CO2 molecules released by shocks from 
the grain surface are likely photodissociated (Figure [3] 
and 111) . A detailed modeling of shocks along the outflow 
walls including the strong protostellar FUV irradiation 
is warranted. 

Of the molecules deemed insensitive to geometry in 
Section mm (H2, H+, NH, CO, H3O+, HCO+ and O), 
the presence of a disk (Figure [5]) show significant differ- 
ences for H+, NH, H3O+, HCO+. This is again due to 
the combination of a lower gas temperature due to the 
higher density in the disk and subsequently faster recom- 
bination rate. The total amount of these molecules thus 
remains approximately the same as in the model Stan- 
dard. 



Figure |6] shows the abundance of the FUV enhanced 
CH+ and the FUV destroyed SH for models Standard, 
b2.5, bl.5 and a20. In the midplane, the models bl.5 
and a20 have a higher abundance of SH by two orders of 
magnitude, while the outflow wall with enhanced CH+ 
abundance is much thinner. The region without FUV 
irradiation but temperature above 100 K is shown in the 
inset of the upper panel. This region coincides with the 
region having an SH fractional abundance larger than 
10"^*^, shown as inset in the lower panel. The region 
with fractional abundance larger than 10^* shows a layer 
structure with the molecule being more photodissociated 
closer to the outflow wall. At larger distances (z = 2000 
and z = 10000), the differences between the models are 
smaller, as found in the previous section for the volume 
averaged abundances. 

The abundances of SH, CH+ and H2O in the Disk 
model are given in Figure [3 The layer with high temper- 
ature and enhanced CH+ abundance (disk atmosphere) 
is thinner than in the outflow wall models in Figure [5] 
because of the higher density in the disk surface. Due to 
the larger density dependence of the SH abundance com- 
pared to H2O and CO2, SH does not reach a fractional 
abundance of lO^'* as in the models 1.5 and a20. The 
water abundance on the other hand is even higher in the 
FUV shielded midplane of the Disk model compared to 
the outflow wall models. This is also seen in Tables H] 
and [5l where the volume averaged abundances of CO2 
and H2O are much more enhanced compared to SH. 

4.3. X-ray driven chemistry 

X-rays affect only the innermost region due to geo- 
metrical dilution oc r~^ (Section [3]). In Table [3] we thus 
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TABLE 5 
Ratios of the volume averaged abundances (Eq. [2J to the model Standard. A chemical age of 5 x 10* yrs and a 



PROTOSTELLAR X-RAY LUMINOSITY OF 10'^ ERG S" 



'^ ARE ASSUMED. THE MASS OVER WHICH IS AVERAGED (M400O AND Af2000o) IS GIVEN 
RELATIVE TO MODEL Standard. 





NH+ 


C+ 


C 


CO 


C02 


H2O 


H2O+ 


H3O+ 


HCO+ 





Within 


4000 AU 


MioQo 


+" 


+ 


(+) 





- 


- 


+ 











Spherical 


2.1 


6(-5) 


3(-5) 


2(-2) 


1.3 


2(1) 


6.8 


l(-4) 


1.6 


5.9 


0.4 


Standard 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Disk 




9.1 


8(-2) 


0.1 


0.2 


1.3 


3(2) 


7(1) 


8(-2) 


9(-2) 


6(-2) 


0.5 


a20 




1.8 


0.3 


0.2 


0.4 


1.3 


4.7 


4.7 


0.2 


0.4 


1.5 


0.6 


a40 




1.5 


0.7 


0.5 


0.8 


1.2 


3.0 


2.5 


0.4 


0.6 


1.5 


0.7 


a 80 




0.64 


1.2 


1.6 


1.1 


0.8 


4(-2) 


0.1 


1.6 


1.7 


0.5 


1.3 


alOO 




0.36 


1.6 


2.4 


1.3 


0.5 


2 (-4) 


0.2 


2.4 


2.8 


0.5 


1.8 


7O 






0.7 


1.0 


1.0 


1.0 


1.2 


1.2 


0.9 


1.0 


1.0 


1.0 


■yO.Ol 






0.6 


1.0 


1.0 


1.0 


1.3 


1.1 


0.9 


0.9 


1.1 


1.0 


-1O.I 






0.1 


3(-2) 


0.9 


1.3 


3.2 


2.5 


2(-2) 


0.4 


2.0 


0.5 


hi. 5 




1.3 


0.6 


0.6 


0.7 


1.1 


2.5 


2.3 


0.6 


0.6 


1.3 


0.8 


b2.5 




0.76 


1.3 


1.4 


1.2 


0.9 


0.2 


0.3 


1.6 


1.4 


0.7 


1.2 


L2e3 






0.2 


0.7 


0.9 


1.1 


0.7 


1.6 


0.6 


1.0 


1.1 


0.8 


Lle4 






0.9 


1.0 


0.9 


1.0 


1.0 


1.3 


1.0 


0.8 


1.1 


0.9 


L4e4 






0.8 


1.1 


1.2 


1.0 


1.4 


1.0 


1.2 


1.1 


0.9 


1.1 


Within 


20000 AU 


Af20000 






















Spherical 


1.8 


l(-3) 


7(-5) 


9(-2) 


1.0 


6.3 


3.5 


l(-3) 


2.4 


2.2 


0.8 


Standard 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Disk 




1.2 


0.5 


0.8 


0.9 


1.0 


5(1) 


6.1 


0.7 


0.7 


0.7 


1.0 


a20 




1.2 


0.4 


0.2 


0.4 


1.0 


3.3 


1.8 


0.2 


1.7 


1.8 


0.8 


■yO.l 






0.1 


4(-3) 


0.1 


1.0 


1.4 


1.3 


7(-3) 


0.9 


1.2 


0.8 


L2e3 






0.2 


0.5 


1.1 


1.0 


1.6 


1.4 


0.2 


0.9 


1.0 


0.9 



Note: a{b) = a X W'' . 

"Ratio to the spherical model within 4000 AU (Table [S] + more than 10-*, 



less than 0.1, (+) more than 3, within 0.1 to 3). 



give the volume averaged abundance within 4000 AU. In 
the two-dimensional model (Standard), the largest dif- 
ferences between the model for an X-ray luminosity of 
10^^ erg s~^ and no X-rays are found for H;^ with an 
X-ray enhancement of a factor of 6. The abundances of 
NH, C and HCO''" are enhanced by a factor of about 
three. Compared to the spherically symmetric model, C 
and HCO+ are enhanced on a similar level, while NH 
is much less enhanced. Destruction of CO2 and wa- 
ter by X-rays has b een found in the sph erical models of 
iStauber et all (p0053 and .Stauber et al.l ()2006. ). respec- 
tively. Since the abundances of H2O and CO2 are already 
decreased by FUV radiation, this effect is not seen in the 
two-dimensional model. 

The extra enhancement of diatomic hydrides due to 
X-rays over FUV processing is much smaller in the two- 
dimensional models compared to the spherical model. As 
discussed in Section 14.1. 1[ the enhancement in the out- 
flow walls by FUV is much larger than by X-rays. Only 
the volume averaged abundance NH, which is X-ray sen- 
sitive but not much enhanced by FUV radiation shows 
some dependence on X-ray irradiation (Figure [3]). How- 
ever note that NH has high energy levels and thus is more 
excited in the outflow walls so that the effect of X-rays 
on the line flux is probably very small. However, a de- 
tailed calculation of the line flux should be carried out 
once collision rates of NH have become available. 

The X-ray enhancement/destruction in different mod- 
els is studied in Table H] by the ratio of the volume av- 
eraged abundance including protostellar X-ray emission 
{Lx = 10^2 erg s'^) and no X-rays within 4000 AU. The 
X-ray enhanced species H^ , NH and HCO+ are similarly 



enhanced for the 2D models except for models a80, a 100 
and Disk where the enhancement in smaller. Models a80 
and a 100 have less mass with only X-rays but no FUV 
irradiation (Section |31) and the molecules are photodisso- 
ciated in the region with X-ray irradiation. In the Disk 
model, however, these molecules are less FUV dissoci- 
ated due to the higher density and even enhanced in the 
warm disk atmosphere. This enhanced layer dominates 
over the X-ray enhanced region and thus the X-ray en- 
hancement is smaller. 

X-r ay destruction of H2O and CO2 f see IStauber et al.1 
120051 [2OO6) is remarkable in models which have a larger 
region with temperature above 100 K but no FUV ir- 
radiation (Models Spherical, a20, a40, jO.l and L4e4)- 
Models with a smaller or similar amount of H2O or CO2 
behave similarly to the model Standard, discussed before. 

We conclude that the influence of X-rays on the chem- 
ical composition is relatively small in this scenario of 
FUV irradiated outflow walls. Mole cules predicted to 
be X-ray tracers (jStauber et al.l [20051 ) and in particular 
diatomic hydrides are much more enhanced by FUV ra- 
diation compared to X-rays. Exceptions are the water 
and CO2 destruction by X-rays which is seen in some of 
the models that are similar to spherical geometry (e.g. 
small outflow opening angle). Also, HCO"'" and H^ are 
considerably enhanced for some models. The differences 
between the models and the uncertainties of the chemi- 
cal model are however as large as the extra enhancement 
by X-rays. For example the total amount of HCO"'" (ob- 
tained from Tables IH [5] and [6]) within 4000 AU is larger 
for the model a40 without X-rays compared to the model 
a80 with X-rays. Thus, it is difficult to establish X-ray 
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Fig. 5. — Density structure and regions witfi T > 100 K and/or FUV irradiation of tfic model Disk. Abundances of CH+, SH and H2O 
are given along cuts of constant z (in the midplane and at 2 = 200, 1000 AU). Tlie model with X-rays (dashed line) assumes an X-ray 
luminosity of lO'''^ erg &~^ . 



TABLE 6 

Enhancement or destruction by X-rays for selected 

SPECIES. The ratio of the volume averaged abundance 

WITHIN 4000 AU between a model with Lx = lO'^-^ erg s~-1 

AND no PROTOSTELLAR X-RAY EMISSION IS SHOWN 



(W4 



= 10^ 



/(^)4000 )■ 



Model 


M4000 


H^ 


NH 


CO2 


H2O 


HCO+ 


Spherical 


2.1 


6.6 


12 


0.08 


0.1 


4.0 


Standard 


1 


5.6 


1.8 


1.0 


1.6 


3.5 


Disk 


9.1 


1.8 


1.1 


0.7 


0.7 


1.1 


a20 


1.8 


5.9 


2.8 


0.01 


0.1 


4.1 


a40 


1.5 


6.0 


2.2 


0.03 


0.4 


4.4 


a 80 


0.64 


4.2 


1.2 


1.5 


1.3 


1.7 


alOO 


0.36 


2.6 


1.0 


1.9 


1.0 


1.0 


7O 


1 


5.9 


2.0 


0.9 


1.5 


3.7 


-yO.Ol 


1 


5.5 


1.9 


0.5 


1.5 


3.7 


-yO.l 


1 


6.6 


1.9 


0.05 


1.0 


5.2 


hi. 5 


1.3 


5.3 


1.6 


0.03 


0.4 


4.0 


b2.5 


0.76 


5.4 


1.6 


1.4 


2.2 


2.6 


L2e3 


1 


6.7 


2.1 


2.1 


1.6 


3.5 


Lle4 


1 


5.6 


1.6 


1.3 


1.6 


3.6 


L4e4 


1 


5.3 


2.0 


0.05 


1.3 


3.3 



tracers for these models with very strong FUV irradia- 
tion. 

In envelopes of low-mass YSOs, however, FUV radia- 
tion could be less dominant compared to X-rays. Several 
effects control the relative importance of FUV and X- 
rays. The colder photosphere of low-mass YSOs emits 
less FUV. Bow shocks, internal jet shocks and accretion 
disk boundary lay er can however still produce a consid- 
erable FUV field (jvan Kcmpcn et al.. i2009a. c). On the 



other hand, most commonly observed low-mass YSOs 
are closer. Since the X-ray luminosity of low-mass YSOs 
is of similar strengths compared to high-mass YSOs (e.g. 
INakaiima et al.|[2003D the X-ray irradiated gas may fill 
the beam and thus the influence of X-rays on observable 
molecular lines becomes stronger. A detailed discussion 
is warranted. Very high resolution interferometric obser- 
vations of FUV and X-ray tracer with ALMA will also 
help to disentangle the effects of FUV and X-ray irradia- 
tion. In particular the X-ray tracer NH, discussed before, 
has a line at 946 GHz within ALMA band 10. High-J 
HCO+ lines (e.g. the J = 10 ^- 9 line at 892 GHz) which 
probe the innermost part of the envelope are also within 
ALMA band 10. 

5. MOLECULAR EXCITATION AND LINE FLUXES 

Observable quantities of the models presented here are 
molecular or atomic line fluxes. In this section, we study 
different excitation effects of diatomic hydrides, like colli- 
sional or radiative pumping. We use a multi-zone escape 
probability method described in Appendix [B] Herschel 
observations are then simulated from the excitation cal- 
culation. For simplicity, we restrict the discussion to 
two prototypical diatomic hydrides, one found to be de- 
stroyed (SH) and enhanced (CH+) in the outflow wall. 
Line fluxes of other species will be presented in an further 
paper together with observations. 

Molecular data for the line frequency and Einstein 
coeffi cients are obtained from the CDMS (iMiiller et al 



2001") for CH+ and the JPL database (jPickett et al 
1998. 1 for SH. The molecular collision rates of many di 
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Fig. 6. — Abundances of CH+ and SH for the models Standard, b2.5, bl.5 and a20. The abundances are given for cuts of constant z 
assuming no protostellar X-ray radiation (solid line) and an X-ray luminosity of 10"^^ erg s~i (dashed line). As an inset to the plots for 
CH+, the regions with FUV irradiation and/or temperature > 100 K are given. The inset to the plots for SH show the abundance of SH 
divided into regions with fractional abundances > 10~*, < 10"^'' and in between. 



atomic hydrides are poorly known and extrapolation is 
th us necessary. For SH-H2, we adopt the OH-H2 rates 
bv lOffer fc van Dishoedd ()1992[ ). scaled for the reduced 
mass. Since collision rates for ortho and para H2 exist, 
we use a thermal ortho to para ratio'*. For CH+-H2 and 
CH+- H, w e scale the CH+-He rat es bv iHammami et aLl 
()2008[) and lHammami et al.l ()2009[ ). As CH+ is abundant 
in regions with electron fraction of or der 10~^, we als o 
consider excitation by electron impact ()Lim et al.lll999[ ). 
A general discussion on the un certainty of collision rates 
is given in iSchoier et al.l (|2005D . 

Transitions observable by Herschel (HIFI/PACS) are 
summarized in Tables [7] and El They give the transition 
wavelength or frequency {vul^ ^ui), the Einstein- A coeffi- 
cient Aui , the energy of the upper level E^i . Transitions 
of CH+ with J^p = 2 - 6 can be observed by PACS, 
while the J = 2 ^ 1 and J = 1 -^ lines at 835.079 
GHz and 1669.170 GHz, respectively, are accessible to 
HIFI. The SH lines are in HIFI bands, except for the 
lines at 1383 GHz, which are given for completeness. For 
SH, only the strongest hyperfine components are used 
and overlapping components have been combined. The 
critical density (ricrit = A^i/ X); ^^1) gives the necessary 
density of collision partners for substantial population of 
the upper level. 

We use wavelength (/xm) for CH+ as most lines are ob- 
servable with PACS and frequency (GHz) for SH (HIFI). 
Levels of CH+ are labeled by the rotational quantum 
number J . The molecular structure of SH is similar to 
hydroxyl (OH) with a ^H electronic ground-state, A-type 
doubling and hyperfine splitting. Levels are thus denoted 
by ^Ho, J, F + /— , with the total angular momentum Q,, 
the rotational angular momentum J, the nuclear spin F 

4 0H2/PH2 = min(3,9 X exp(-170.6/Tgas)) 



TABLE 7 
Molecular data of SH. The critical density is given for a 

TEMPERATURE OF 100 K AND A THERMAL ORTHO TO PARA RATIO 
OF H2. 





l^ul 


Aui 


Eui 


ncrit(H2) 


[J,F/Parity] 


[GHz] 


[s-i] 


fKl 


[cm-3] 


2Hi/2 3/2, 1- ^ 1/2, 1-f 


866.947 


1.9(-3) 


571 


4.3(6) 


2Hi/2 3/2,2+^1/2,1- 


875.267 


1.5(-3) 


572 


4.7(6) 


2H3/2 5/2, 3+ -i- 3/2, 2- 


1382.910 


9.4(-3) 


66 


3.8(7) 


2H3/2 5/2,3-^3/2,2-1- 


1383.241 


9.4(-3) 


66 


4.7(7) 


2Hi/2 5/2,2+^3/2,1- 


1447.012 


1.6(-3) 


641 


2.6(7) 


2Hi/2 5/2,3- -s-3/2,2+ 


1455.100 


1.6(-3) 


642 


3.0(7) 


2H3/2 7/2,4-^5/2,3+ 


1935.206 


3.5(-2) 


159 


1.2(8) 


2H3/2 7/2,4+^5/2,3- 


1935.847 


3.5(-2) 


159 


1.4(8) 



TABLE 8 

Molecular data of CH+. The critical density is given for 

A temperature of 500 K. 





^ul 


Aui 


Eui 


ncrit(H2) 


«crit(H) 


"crit(e ) 


fJl 


[/im] 


fs-M 


fKl 


[cm-3] 


[cm-3] 


[cm-3] 


1 ^0 


359.0 


6.4(-3) 


40 


5.3(7) 


3.8(7) 


3.0(4) 


2-> 1 


179.6 


6.1(-2) 


120 


2.5(8) 


1.9(8) 


1.3(5) 


3->2 


119.9 


2.2(-l) 


240 


7.1(8) 


5.2(8) 


3.8(5) 


4-^3 


90.0 


5.4(-l) 


400 


1.6(9) 


1.2(9) 


7.8(5) 


5^4 


72.1 


1.1(0) 


600 


3.2(9) 


2.3(9) 


1.7(6) 


6^5 


60.2 


1.9(0) 


838 


5.8(9) 


4.2(9) 


3.0(6) 



and the parity (+/-). The CH+ and SH molecules show 
typical properties of diatomic hydrides, with high criti- 
cal densities (> 10® cm"'^ for H2 collisions) and large en- 
ergy separation between subsequent levels. This is due to 
the large rotation constant combined with a large dipole 
moment. Line frequencies thus mostly lie above the at- 
mospheric windows. Excitation mechanism other than 
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coUisional excitation are likely. For example, pumping 
by dust continuum radiation or formation in an excited 
state have been suggested by iBlackl ()1998D for CH+ . 

5.1. Excitation: Level population 

The calculated level populations of SH and CH+ at 
different positions of model Standard are given in Fig- 
ure [71 For SH, the position in the midplane with peak 
abundance and at the edge of the cloud are presented. 
For CH"^ , we show positions with peak abundances along 
cuts through the midplane and at z = 2000 AU and 
z — 10000 AU (Section l4J|) . In order to study the impor- 
tance of different excitation mechanism, we give the nor- 
malized level populations depending on the upper level 
energy for the following models: 

• All. CoUisional and radiative excitation considered. 

• No Dust. Dust radiation switched off by setting the 
dust opacity to zero. This model allows to study 
the importance of dust pumping. 

• No Radiation. Ambient radiation field switched off 
{{Jij) = 0). This model together with the previous 
allows to study the importance of pumping by line 
emission. 

• No Collisions. Density of the collision partners set 
to zero. This situation corresponds to pumping 
purely by dust radiation. 

• LTE. Level population in the local thermal equilib- 
rium (LTE) with the gas temperature. This model 
is given for comparison as the assumption of LTE 
is often made due to missing collision rates. 

• Tform = 100 K, Tfoim = 3000 K (Only CH+). Col- 
lisional excitation, radiative excitation and excited 
formation considered with formation temperature 
Tform (Appendix [B| . We adopt a Tfomi of order 
of the kinetic temperature (3000 K) or much lower 
(100 K). 

For SH, the population of the hyperfine quadruplet is 
combined. We note that the level population can only 
be translated directly to emissivity per volume in the case 
of optically thin radiation. Figure [8] gives different com- 
ponents (line emission, dust emission, cosmic microwave 
background and stellar photons) of the ambient radiation 
field (Jy) at positions A (SH) and C (CH+). 

The CH+ level population is far out of LTE. It peaks in 
the J = 1 state in position C and the ground state in po- 
sitions D and E. For all positions, the population rapidly 
decreases to higher J levels and from position C to E for 
the same level. CoUisional excitation (model No Colli- 
sions) governs the excitation with increasing differences 
to model All for higher J . A model without electron 
excitation but collisions with H and H2 has level popula- 
tions within 20 % to the model All. Electron excitation 
is thus unimportant relative to other effects. Models No 
Dust and No Radiation differ less relative to model All in 
CH"*" compared to SH, thus radiation is less important. 
Unlike SH, the models No Dust and No Radiation differ, 
due to the contribution of line photons to the ambient 
radiation field. At position C, line photons dominate the 



radiation field up to the J = 5 — > 4 transition (Figure 
[8] right). Optical depths of up to r « 1 (line and dust) 
for vertical cuts through the outflow wall are reached in 
the lower J transition. For levels connected by these 
transitions, we find the largest differences between mod- 
els No Dust and No Radiation. We conclude that CH+ 
is mainly collisionally excited with small contributions 
from pumping by line trapping and dust radiation. 

The importance of collisions for the excitation of CH+ 
despite the density of the collision partners below the 
critical density can be understood by the high kinetic 
temperature. Radiative transitions connect only levels 
with AJ = 1 due to selection rules. Rates for radia- 
tive pumping are proportional to the lower level popula- 
tion of a pumping transition and thus decrease quickly 
to higher J. CoUisional excitation is not restricted to 
AJ = 1 and can pump over several levels, especially 
for high kinetic temperatures, when the upward rates 
are similar or even higher than the downward rates 
{Ciugi = C„/5„exp(-Ai;/A:Tki„)). 

At position C in the model Standard, the radiative ex- 
citation rate decreases from about 10^^ ( J = 1 ^- 2) to 
10-5 cm-3 s-i (J = 9 ^ 10), for Bi^{J^i) and S„i(J„i) 
of order 0.1 s^^ and a CH"*" density of 0.4 cm-"^. The col- 
lisional excitation rates from the most populated levels 
(J < 3) to a higher level is found to be about 5 x 10""^ 
s~^ cm""^ into J = 4 and slightly less for higher J, for 
rate coefficients Kui of order a few times 10^^*^ cm'^ s~^ 
and a density of 10^ cm^^^. Thus, for J — A, radiative 
and coUisional rates are similar, but collisions dominate 
for higher J. This is in agreement with the findings in 
the previous paragraphs. 

Excited formation affects the level population rela- 
tively little with considerable differences only for higher 
J levels at lower densities (position E) and a high forma- 
tion temperature (Tform = 3000 K). At position C, the 
total formation rate of CH+ is 6 x 10"'^ s""'^ cm"*^, larger 
than the coUisional excitation, but the effect on the level 
population is small due to the distribution over differ- 
ent levels following the Boltzmann statistic. Due to this 
distribution, which supports the thermalization to the 
formation temperature, the infiuence of excited forma- 
tion is mostly seen in higher J levels for Tform — 3000 K, 
where coUisional excitation competes less with molecules 
formed in an excited state. 

Observation of the J = 2— >luptoJ = 6^'5of 
CH+ using the In frared Space Obs e rvator y (ISO) have 
been reported by iCernicharo et al.] (J1997I ) towards the 
planetary nebula NGC 7027. They can be explained by 
a single excitation temperature of 150 ±20 K for a kinetic 
temperature of about 300-500 K and densities of a few 
times 10^ cm~^. It is interesting to note that, at posi- 
tion C of our model, the level population up to J = 6 can 
be well fitted by a single Boltzmann dist ribution with a 
temperature of 129 K. Note that Cernicha ro et al.l ()1997f ) 
have used N2H+ collision rates and a slightly higher den- 
sity. Our calculation thus suggests that the observed line 
ratio can be explained even with a much higher kinetic 
temperature as is required for the formation. A more de- 
tailed model of their source is however required to draw 
further conclusions. 

The SH level population is out of LTE like CH+. The 
molecule is however destroyed by FUV radiation and thus 
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Fig. 7. — Normalized level population of SH (top panels) and CH"'' (bottom panels) at different positions. The positions A-E are given 
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Fig. 8. — Radiation field of SH and CH+ at positions A and C, 
respectively. The positions are given in Figure [7] 

abundant in regions with lower temperatures compared 
to CH+ . The level population is thus concentrated in the 
ground state for both positions A and B. Levels of the 
^ni/2 ladder with level energy above 500 K are sparsely 
populated with relative population below 10^"*. The 
agreement of the model No Collisions with the model All 
in the ^113/2 states indicates that these level are pumped 
by dust radiation. The '^Ili/2 states on the other hand 
are connected by lines with lower frequency and thus 
weaker dust radiation field. These levels are pumped by 
both collisions and dust radiation. 

5.2. Line fluxes 



Calculated line fluxes of SH and CH+ for different mod- 
els are given in Table El Synthetic frequency integrated 
maps and spectra at selected positions for the J = 1 -^ 
and J = 6 — ;■ 5 lines of the model Standard are shown in 
Figure ini Maps and spectra are obtained from the level 
population discussed in the previous section. Synthetic 
maps are calculated by integrating the radiative transfer 
equation. To simulate observations, the maps are con- 
volved with the appropriate Herschel beam, centered on 
the source and assumed to be a Gaussian. For the ray- 
tracing, we assume the same physical parameters as in 
paper II. These are a distance of 1 kpc, an inclination of 
30° and a position angle of the outflow direction of 96° . 

The predicted CH+ fluxes are remarkably strong, even 
for higher J lines. The integrated intensities correspond 
to fluxes of order a few times 10^^^ W cm~^ and are 
easily detectable with PACS or HIFI. Note that they 
are similar to those reported towards NGC 7027. Strong 
lines for higher J transitions despite the low level popula- 
tion found in the previous section are due to the smaller 
beam at shorter wavelengths and the emissivity (in K 
km s^^) per excited molecule oc Aui/v'^i, increasing by a 
factor of 10 from the J = 1 ^> to the J = 6 ^^ 5 tran- 
sition. Using a Boltzmann diagram, corrected for line 
absorption, the model Standard fits well for a column 
density of 7 x lO""^^ cm~^ and an excitation temperature 
of about 148 K {Slab at LTE). 
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TABLE 9 
Molecular line fluxes of SH and CH+. 



Model 




Line Flux 
















[K km s-i]* 


[K km s-i] 


[K km s-i] 


[K km s-i] 


[K km s-i] 


[K km s-i] 


CH+ 




1 -^0 


2^1 


3-^2 


4^3 


5^4 


6^5 


Wavelength [/jm] 




359.0 


179.6 


119.9 


90.0 


72.1 


60.2 


Herschel-Beam ["] 




25.8" 


12.9" 


8.6" 


6.5" 


5.2" 


4.3" 


Standard"' 




12.3 


23.9 


20.2 


10.5 


4.6 


2.1 


Slab at LTE^ 




9.8 


22.5 


22.4 


13.6 


5.3 


1.6 


Standard (2 kpc) 




4.3 


8.9 


8.0 


4.2 


1.9 


0.9 


Standard (LTE) 




12.0 


113.1 


333.0 


612.3 


871.9 


1064.9 


Standard (Ttar-m = 


100 K) 


14.9 


26.4 


21.8 


11.0 


4.7 


2.1 


Standard (Ttarui = 


3000 K) 


15.8 


30.8 


29.4 


18.1 


9.7 


5.4 


Standard ((P ) 




14.8 


29.5 


26.1 


14.1 


6.3 


3.0 


Standard (90°) 




16.1 


23.6 


17.3 


8.0 


3.2 


1.3 


Standard (25.8") 




12.3 


7.5 


2.9 


0.9 


0.3 


0.1 


Standard (Noturb) 




9.4 


19.0 


19.3 


11.2 


5.0 


2.2 


Disk 




8.2 


11.7 


7.3 


2.6 


0.8 


0.3 


a 20 




4.9 


10.1 


10.7 


7.2 


3.5 


1.6 


alOO 




28.8 


40.7 


29.8 


14.6 


6.7 


3.3 


SH^ 




^ni/2(3/2^1/2) 


^n3/2(5/2^3/2) 


ni/2(5/2-»3/2) 


^ni/2(5/2-»3/2) 


^n3/2(7/2^5/2) 


n3/2(7/2->5/2) 


Frequency [GHz] 




866.947 


1382.910 


1447.012 


1455.100 


1935.206 


1935.847 


Herschel-Beam ["] 




24.9" 


15.6" 


14.9" 


14.8" 


11.1" 


11.1" 


Standard'^ 




3.8(-5) 


-2.1(-3)'' 


4.6(-5) 


5.0(-5) 


5.8(-4) 


4.4(-4) 


Standard (2 kpc) 




3 


2(-6) 


-9.7(-4) 


2.1(-6) 


2.4(-6) 


8.7(-5) 


6.1(-5) 


Standard (LTE) 




8 


2(-5) 


3.2(-l) 


8.8(-4) 


8.5(-4) 


9.1 (-2) 


9.0(-2) 


Standard (OP) 




8 


4(-5) 


5.4(-3) 


1.2(-4) 


1.3{-4) 


9.2(-4) 


7.3(-4) 


Standard (90P ) 




9 


9(-6) 


-6.1(-2) 


5.7(-6) 


6.6(-6) 


-7.9(-4) 


-7.6(-4) 


Disk 




1 


8(-5) 


-9.9(-4) 


2.4(-5) 


2.6(-5) 


2.8(-3) 


2.3(-3) 


a 20 




1 


6(-4) 


-1.5(-1) 


1.8(-4) 


2.0(-4) 


3.8(-2) 


3.0(-2) 


alOO 




8 


2(-5) 


1.7(-3) 


9.8(-5) 


l.l(-4) 


2.8(-4) 


2.1(-4) 



Note: a{b) = a X lO'' . 

"If no other indication is given, the inclination is 30° and the distance 1 kpc (see Appendix [Bj 

''The conversion factor from K km s~^ to W cm~^ is 2.7 X 10~^'^ / (\[fj,m])^ . 

'^Fit to model Standard for a CH+ column density of 7 X 10^^ cm~^ and Tex = 148 K. 

''Negative fluxes mean line in absorption. 

'^See Table [7] for the exact label of the transition including fine structure. 



Different radiative transfer models of CH+ derived 
from the chemical model Standard have similar fluxes in 
the J = 1 — >■ line, with deviation smaller than about 30 
%. An exception is the model Standard (2 kpc) which as- 
sumes a distance of 2 kpc instead of 1 k pc to the source, 
as suggested by [Schneider et al.l (J2006D . For higher J 
transitions, fluxes however differ significantly between 
the models. Assuming the level population in the LTE 
[Standard (LTE)) yields fluxes more than two orders of 
magnitude larger compared to the model Standard. Ex- 
cited formation increases the line fluxes of all lines by less 
than 20 % for Tform = 100 K. Assuming Tform = 3000 K, 
a larger increase is found for higher J lines and amounts 
to a factor of 3 for the J = 6 — )> 5 line. Inclination affects 
the line flux by the combination of line and dust opacity, 
together with the beam containing different parts of the 
model. For example, the dust opacity reaches about 1.2 
for the J = 6 — )► 5 line, while it is below 0.07 for the 
J = 1 ^> line. Thus, the flux of the J = 6 — >■ 5 line in 
the model Standard (9(F) (edge on) is about a factor of 
two smaller compared to the model Standard ((F) with 
a line of sight parallel to the outflow. Compared to the 
model Standard, the inclination changes the flux mostly 
by less than 40 %. We note that for other sources with 
larger total column density, the dust attenuation may be 
important even for longer wavelengths. Increasing the 



distance to 2 kpc flnally, decreases the fluxes by a factor 
of about 2.5. 

The frequency integrated maps and spectra of model 
Standard given in FigurelHlshow that the concentration of 
the emission to the innermost region is more pronounced 
for the J = 6 —> 5 line compared to the J — 1 ^ 0. This 
is also reflected in a model (Standard (25.9")), convolved 
to a 25.8" beam. The flux ratios of that model to model 
Standard decreases almost as fast as the beam dilution 
factor, indicating that most of the emission is within the 
beam. The profile of the J = 1 — > line shows con- 
siderable self- absorption in the region where the line of 
sight follows the outflow wall due to the opening angle 
being similar to the inclination (positions A and B). The 
optical depth at the line center at these positions are 7 
(A) and 19 (B). Further out, the optical depth is about 
3 (C), which is about twice the value found in the Sec- 
tion 15.11 for a ray crossing the outflow wall once. The 
optical depth is even small for distances larger than po- 
sition C. The J = 6 —^ 5 line has much smaller line 
optical depths below 0.2 (position A) compared to the 
J = 1 — i> line. Beam convolved spectra of both lines 
however, do not show a self-absorption dip, suggesting 
that most of the emission comes from regions with mod- 
erate optical depth. 

The considerable line optical depths for lines of sight 
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Fig. 9. — Synthetic velocity integrated maps (top panels) and spectra (bottom panels) of the CH+ J = 1 — > (left panels) and J = 6 — ^ 
transitions (right panels) for the model Standard. The appropriate Hcrschel beam (FWHM) is given by a circle. Contour lines give 
integrated intensities of 10 and 100 K km s~ . The continuum is subtracted from the spectra and the red (gray) spectrum is convolved to 
the Herschel beam and multiplied by 10. The spectra at positions A, B and C are given in a pencil beam. 



along the outflow walls raises the question, how much 
the modeled line fluxes depend on the assumed intrin- 
sic line shape and velocity structure. HIFI can observe 
the J = 1 ^ and J = 2 -> 1 hues of CH+ spectrally 
resolved but with a beam corresponding to more than 
10000 AU at the assumed distance of 1 kpc. This is 
insufhcient to directly probe the velocity structure along 
the outflow wall. It can however be used to constrain the 
average line width along the outflow wall. Spatially re- 
solved interferometric observations reveal broader lines 
in the outflow wall wi thout much veloc ity shift to the 
systemic velocity (.Bruderer et al.l[2"009al ). This has been 
implemented in model Standard ( Appendix [B|) . A model 
without increased microturbulence in the outflow walls 
{Standard (Noturb)) has been calculated and shows to- 
tal line fluxes within 20 % to model Standard. We con- 
clude that the while the modeled line shape may depend 
strongly on the assumed intrinsic line and velocity struc- 
ture, the line fluxes are pretty independent. 

Different chemical models {a20, alOO and Disk) affect 
the line fluxes of CH+ by a combination of the chemical 
abundance, the temperature and density structure and 
also the beam containing different regions of the model. 
For example the total amount of CH+ in the models a20, 
a 100 and Disk within 4000 AU relative to the model 
Standard are 0.64, 0.36 and 1.04, respectively (Tables |4] 
and [5]) . The line fluxes relative to model Standard are 
thus similar to the abundance ratios for models a20 and 
Disk in low- J lines. For a 100, the line is however much 
stronger due to the larger mass with high temperature 



(Table [2]). For higher J lines, the situation is different. 
The lines of the model Disk get weaker with J relative to 
model Standard, while those of model a20 get stronger. 
This is due to the fact, that the abundance of CH+ de- 
creases in model Disk in regions with high density, while 
in the model a20 more material with high density is 
heated by FUV. We conclude that the low- J lines of CH"*" 
are better tracers for the chemical abundance than higher 
J lines, as the higher J lines depend more on the tem- 
perature structure. The difference between the models 
is of order of about 3, similar to the uncertainty of the 
chemical model. Details of the geometry can thus not be 
obtained solely by observations of diatomic hydrides. 

The SH lines are very weak, with most line fluxes 
clearly below the detection limit. This is remarkable 
since the volume averaged abundance of SH (Table [3]) 
is only about two orders of magnitude smaller than for 
CH+. The weak fluxes can be explained by the poor ex- 
citation in regions with high abundance of SH, due to 
the large energy separation of the levels (Section 15. ip . 



Lines of the II1/2 ladder with upper level energies of 
several 100 K have line fluxes below 10""^ K km s~^. 
Lines between ^113/2 states can reach line fluxes of a few 
times 10~^ K km s"-'^ owing to their level energies simi- 
lar to the dust or gas temperature in regions where SH is 
abundant. These lines can also be seen in absorption if 
a sufficient column density of cold dust lies between the 
SH emitting region and the observer. Unfortunately, this 



feature is most prominent in the H; 



3/2(5/2^3/2) 



line at 
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1383 GHz which cannot be observed by PACS or HIFI. 
The abundance of SH within 4000 AU of the models 
a20, alOO and Disk relative to model Standard are 180, 
0.04 and 23 (Tables H and EJ. This is only reflected in 
lines of the ^113/2 levels which are excited in larger parts 
of the envelope compared to the ^lli/2 levels. Because of 
the larger amount of warm material, lines obtained from 
model a 100 can be stronger than from the model Stan- 
dard despite the lower abundance. We conclude that the 
excitation in the warm region of the FUV outflow walls is 
essential for observable lines of diatomic hydrides. Thus, 
observations or non-detections do not necessarily reflect 
column densities, due to high critical densities combined 
with large energy separations between subsequent energy 
levels. 

6. OBSERVING DIATOMIC HYDRIDES WITH 
HERSCHEL AND ALMA 

The previous sections have shown that light diatomic 
hydrides (CH, NH, SH, OH) and their ions belong to a 
particular class of molecules that are expected to probe 
the outflow walls of YSO envelopes. This is due to three 
factors relating to formation and excitation of the di- 
atomic hydrides: i.) Owing to considerable activation 
barriers for the formation of many of the diatomic hy- 
drides, a high temperature of several 100 K is needed for 
their formation, ii.) The chemical network leading to 
the formation is initiated by photoionization and large 
FUV fields are required for the formation of diatomic 
hydrides, iii.) Due to high critical densities and a large 
energy separation between the molecular levels, a high 
density and temperature is needed for excitation of the 
diatomic hydrides. 

Protostellar FUV radiation can escape through a low 
density cavity, etched out by the outflow, and irradiate 
the outflow walls separating the outflow and envelope. 
Our work predicts the abundance and excitation of di- 
atomic hydrides in this thin outflow walls to govern the 
total line flux and result in fluxes that can be easily de- 
tected. Diatomic hydrides thus are predicted to be valu- 
able tracers for the chemistry and physics of the outflow 
walls, which are in many ways similar to the upper at- 
mosphere of protoplanetary disks. 

The Herschel Space Observatory provides a unique op- 
portunity to observe diatomic hydrides, as their lines are 
above atmospheric windows. Beyond the relatively large 
sample of hydride lines in the HIFI spectral range, which 
is targeted by the WISH Herschel guaranteed time key 
program in this source, PACS observations of higher J 
lines may yield important constraints on the hottest re- 
gions. A particular example are high-J lines of CH+ 
which are predicted to be detectable by PACS. For this 
particular source, the relatively large beam of Herschel 
does not allow to spatially resolve the outflow walls and 
our model do also not indicate that the predicted line 
fluxes can be used easily as tracer for the exact shape of 
the cavity. In other sources with a good balance between 
distance and affected outflow area (e.g Serpens-FIRS 1), 
the good spatial resolution of the PACS camera may al- 
low to spatially resolve different transitions of CH+ . 

The upcoming Atacama Large Millimeter and submil- 
limeter Array (ALMA) can observe diatomic hydrides 
at very high angular resolution using band 10 receivers 



(787-950 GHz). For example the J = 1 -> line of the 
^•^CH"'' isotope is s hifted sufficiently to b e inside an at- 
mospheric window (Falgarone et al.ll2005| ) . Also 0H+ at 
909 GHz (jWv rowski ct al. 2010) and SH+ at 893 GHz 
(jMenten et al.l I201Q) have been detected from ground. 
Furthermore, lines of NH (946 GHz) and SH (866 and 
875 GHz) are observable with ALMA band 10 receivers. 
Herschel observations can thus be complemented in the 
near future with high-resolution interferometric data to 
constrain parameters like the geometry in the innermost 
region. 

7. LIMITATIONS OF THE MODEL 

The scenario of directly irradiated outflow walls by ra- 
diation of the protostar is not the only possible to al- 
ter the chemical composition in the outflow walls. For 
example mixing between the warm and ionized outflow 
with t he envelop e has been suggested by iTavlor fc Ragal 
(fT995[) (see also iRaga et"aLl [20021 [2001 . The observed 
abundance of CO"*", however, cannot be explained by 
mixing solely (paper II). Major dissociative shocks may 
occur at the boundary layer between outflow and inflow 
and produce FUV radiation that irradiates the quiescent 
envelope. To replace the FUV field of direct protostellar 
radiation with shocks, a rather high shock velocity is re- 
quired. For example, in the model Stand ard at z = 2000 
AU, f.s ~ 75 km s^^ is required (Neuf eld fc Dalgarnd 
Il989f ) . Such conditions are clearly met for the bow shock 
or internal jet working surfaces. Along the outflow walls, 
however, the velocity perpendicular to the quiescent en- 
velope is likely smaller. In the simple model of the Cou- 
ette flow solution (e.g. iLizano fc Giovanardilll995[ ) ne- 
glecting turbulence, the mean velocity in the entrain- 
ment/mixing is along the outflow cavity and depends 
linearly on the distance to the cavity edge. The width of 
this layer and the substructure of the velocity including 
turbulences are however not easily constrained. Obser- 
vationally, no clear signs of shocks in the outfiow walls 
have bee n found in inter f eromet ric observations of AFGL 
2591 bv iBruderer et al.l (|2009 al We conclude that fur- 
ther study of the combined effects of mixing and shocks 
is warranted, especially in the context of low-mass star 
formation with less or cooler protostellar FUV radiation. 

Several simplifications have been made to make the 
calculation feasible. As discussed above, the physical 
structure does not account for shocks and mixing in the 
outflow layer or a clumpy surface of the outflow wall. The 
chemical models used here do not consider non-thermal 
desorption (e.g. photodesorption) . The photodesorption 
yields are not well known for many species. However, 
the change in abundance at the 100 K evaporation tem- 
perature of H2O, CO2 and H2S is small in the presence 
of a strong FUV field such as in the outflow wall. Thus 
the photodesorption would not affect the results of the 
outflow wall enhanced species considerably. The drop 
of the water abundance in the region with FUV irra- 
diation but temperatures below 100 K would be miti- 
gated by photodesorption. FUV dissociation and ion- 
ization processes are treated here in the rate approach 
which does not allow to specify the protostellar radiation 
fleld explicitly. The surface temperature of the proto- 
star is however not well constrained by observations and 
would affect the chemistry if the temperature were not 
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hot enough to dissociate CO (jSpaans et aLlll995l) . The 
models do not account for y ibrationally excited H2 (e.g. 
[Sternberg fc Dalgarnoll995l ) which may help to overcome 
the activation barrier and further increase the abundance 
of e.g. CH+ or OH. Grain surface reactions, except for 
H2 formation and charge exchange with grains, are not 
considered in the models. For species with longer chemi- 
cal time scales (chemical clock s, e.g. many sulfur bearing 
species, iWakelam et al.l [20041) the effects of infall must 
be taken into account. 

8. CONCLUSIONS AND OUTLOOK 

We have used a detailed chemical model of the high- 
mass star- forming region AFGL 2591 to study the effect 
of different shapes of a concave low-density cavity on the 
chemistry. The cavity, etched out by the outflow, al- 
lows FUV radiation to escape the innermost region and 
irradiate the walls between outflow and envelope. We 
study especially those molecules that are newly observ- 
able with the HIFI/PACS instrument onboard Herschel 
and in particular light diatomic hydrides (CH, NH, SH, 
OH). The model used in this work is based on the two- 
dimensional ax i symme tric chemical model introduced in 
I Bruderer et al.l ()2009bf ) , extended by a self-consistent cal- 
culation of the dust temperature and an escape proba- 
bility method to calculate the molecular excitation. The 
model considers protostellar FUV and X-ray irradiation 
an d makes use of t he grid of chemical models presented 
in I Bruderer et al.l ()2009c| l . The main conclusions of the 
work are: 

1. Protostellar FUV radiation can enhance the abun- 
dance of light diatomic hydrides and their ions by 
photoionization/dissociation processes and heat- 
ing in the outflow-walls considerably. The overall 
amount of CH+ , 0H+ and NH+ in the envelope is 
orders of magnitude larger than in spherically sym- 
metric geometry (Section [4]). Similarly enhanced 
are 0+ and H2O+. These species thus appear to 
be clear tracers of extended FUV irradiation and 
thus the geometric structure of the inner regions of 
YSOs. 

2. Depending whether the vertical height of the en- 
velope from the midplane to the outflow (geomet- 
rically thin or thick), the size and mass of the 
infrared heated hot-core region with temperature 
above 100 K, but no FUV irradiation, can change 
substantially (Section |3]). Geometrically thin inner 
regions allows FUV radiation to penetrate to the 
midplane and destroy the bulk of important species 
such as CO2 and H2O. The shape of the cavity is 
thus also important for species which are destroyed 
by FUV radiation (Section [3J] and gj]). 

3. Provided that FUV radiation can escape the in- 
nermost region, the effects of different geometries 
(opening angle of the outflow, cavity shape, density 
in the cavity) on the total amount of species en- 
hanced in the outflow- waU (CH, CH+, 0H+, NH+, 
C+ and H2O+) are relatively small. Deviations are 
typically less than a factor of three (Section 14. 2p . 

4. In outflow walls, the influence of protostellar X- 
ray irradiation on the chemistry and in particular 



on the diatomic hydrides is relatively small with 
respect to the much larger enhancement by FUV 
radiation (Section 143)) . 

5. The large separation between the energy levels of 
diatomic hydrides requires a hot and dense gas for 
excitation. Thus, the line fluxes of CH+, predicted 
to be enhanced in the hot outflow wall, reach fluxes 
of several K km s~^ even in the J = 6 — >■ 5 line. 
On the other hand, SH, is predicted to be destroyed 
in the outflow wall, but to be abundant in the en- 
velope. Thus, despite the volume averaged abun- 
dance of SH being only two orders of magnitude 
lower than for CH+ , it is not sufficiently excited in 
the cool envelope to be detectable (Section [5?T] and 

[El. 



6. Excited formation of CH+ does not considerably 
affect observable line fluxes. Exceptions are high- J 
lines, if a large formation temperature is assumed 
and the ambient density is low (Section 15.1! and 

[El. 



This suggests that diatomic hydrides are important 
tracers of warm and FUV irradiated material. Di- 
atomic hydrides are chemically closely related to im- 
portant species like water. Studying diatomic hydrides 
will provide important insight in the chemistry of YSO 
envelopes. Not only the chemistry of the diatomic hy- 
drides is interesting, as they may ultimately also be used 
a tracers for physical conditions (like the FUV radia- 
tion). Herschel will for the first time observe diatomic 
hydrides with a good sensitivity and angular resolution 
and making these tracers available. ALMA and in partic- 
ular band 10 will allow to continue the study of diatomic 
hydrides at very high angular resolution. Since diatomic 
hydrides have not been studied thoroughly in the past, 
much molecular data and in particular collisional excita- 
tion rates are missing. More studies on these processes 
will be necessary in order to analyse Herschel observa- 
tions. ^^^^^^^^^^ ^_^^ 

In this series of papers (jBruderer et al.l [2009cl lbl and 
this work), we have introduced a fast method for chem- 
ical modeling of envelopes of young stellar objects and 
studying the molecular excitation. The numerical speed 
up allows to construct more detailed models including 
multidimensional geometries. The method will also fa- 
cilitate the study of larger samples of observations as 
they will soon be available by new facilities with high 
sensitivity and larger spectral coverage as for example 
the upcoming Atacama Large Millimeter/submillimeter 
Array. Ultimately, extensions of the methods presented 
here may also be useful to study the chemical/physical 
evolution of protoplanetary disks. 
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has simplified the development of the escape probabil- 
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supported by the Research Corporation and the NASA 
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APPENDIX 

A. DUST RADIATIVE TRANSFER 

The dust temperature in steady state condition is ob- 
tained by equating the total dust emission Femit (erg s~^ 
g~-^) with the total absorption Fabs (erg s~-^ g~"^)- Dust 
may absorb either stellar photons or ambient emission by 
dust. Assuming a single dust temperature Toust; we can 
write the energy balance at one position of the envelope 
as 



abs 



/ / Kxh{n)dXdn ^An f KxBxiTdust 



)dX 



with K\ (cm^ g 

mass, I\{^) (erg s~^ cm^^ 



(Al) 
the dust absorption opacity per 
/^m~^ sr~^) the intensity de- 
pending on the direction and wavelength and B\{TDust) 
the Planck function. To solve this coupled problem, 
we implement the Monte Carlo approach presented by 
IBiorkman fc Woo"dl ()2001D . The results of the new code 
have been verified with the benchmark tests suggested by 
llvezic et al.l (|1997[) and lStam atellos fc WhitworthI ()2003f ) 
and also by comparing the results with the D USTY code 
(jlvezic fc Elitzuiill997l: iNenkova et al.lll999( ). 

For t he models of AFGL 259 1 , we a dopt the dust opac- 
ities bv lOssenkopf fc HenninsJ ()1994l ) (column 5). This 
dust opacities have been found to explain the observa- 
tion of AFGL 2591 well and yield a "standard" dust 
to m olecule mass ratio of about 100 (|van der Tak et all 
|1999() . For the calculation of the FUV intensity, the dust 
opacities are extended to shorter wavelengths (A < l^m) 
with t he absorption and scattering properties by IDraind 
()2003() . available on his web-page (similar to paper II). 
An example of calculated dust temperature is given in 
Figure [To] for the model used in paper II. The figure also 
shows the calculation domain which consists of 5 nested 
grids. 
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Fig. 10. — Inner region of the calcalultion grid and the 
temperature of an axisynimetric model of AFGL 2591. 
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The dust radiative transfer code also calculates the 
local mean intensity of scattered and attenuated stel- 
lar photons with a Mon te Carlo method similar to 
Ivan Zadelhoff et"aLl (|2003f) . The analytical expression 
for the photodissociation and ionization rates (e.g. 



Ivan Dishoeckl fl988l ) implemented in paper I require the 
calculation of the unattenuated FUV flux Go [ISRF] and 
the visual extinction Ay- We estimate these quantities 
by the mean intensities at both edges of the FUV band 
(6 and 13.6 eV) and at 1 micron. We have tested the 
approach by comparing photoionization and dissociation 
rates obtained fr om the fit and calcu l ated by the FUV 
cross sections of Ivan Dishoeck et al.l ()2006t )^. We find 
agreement within a factor of two in the chemically rele- 
vant range of attenuation {Ay < 10). This is reasonable 
provided that we do not implement the same dust prop- 
erties and spectral shape of the emitting spectra (see Sec- 
tion [7] for a discussion of the spectral shape). The FUV 
calculation is also tested against the results in paper II. 
We find good agreement, for example the mass and vol- 
ume of material irradiated with an attenuated fiux larger 
than 1 ISRF agree to within 30 %. 

B. MOLECULAR EXCITATION ANALYSIS 

For the calculation of the level population of SH and 
CH"*", a multi-zone escape probability method is used. 
This appendix briefly discusses the method. Further de- 
tails and the code will be published elsewhere. 

The excitation of a molecule or atom (i.e. the popula- 
tion of different energy levels) is controlled by collisional 
and radiative processes. The species can be excited or 
deexcited by collisions with H2, H and electrons or by 
radiation either from the dust continuum or molecular 
emission from the same or other molecules. The popula- 
tion Tii (cm~^) of level i at one position is desc ribed by 
the rate equations (e.g. Ivan der Tak et al]|2007D 



dn 



N 



N 



-1T = J2 "J"^:'^ - "* IZ ^^^ +^'- "^^« = ' (^1) 



dt ^-'^'^^ -^-y 



where N is the number of levels considered, Fi (cm~^ 
s""'^) and Di (s~^) are the chemical formation and de- 
struction rates and Pij (s^^) is given by 



P,. 



A,j + B,,{j,j) + a. 



Bij {J, 
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{E, > Ej) 
{E, < Ej) 
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Here, Aij (s^^) and Bij (erg~^ cm^ Hz) are the Ein- 
stein coefficients for spontaneous and induced emis- 
sion, respectively. The collisional excitation and deex- 
citation rates C^ (s~^) are obtained from the sum of 
ncoiKij{T\iin) for different collision partners, with the 
density of the collision partners ricoi and the collision rate 
coefficients i^ij(Tkin) (cm"^ s^^) at a kinetic temperature 
Tkin- The level population is assum ed t o be constant 
with time and thus drii/dt = in Eq. IBll 

The level distribution at formation may govern the 
level population if a molecule is destroyed in collisions 
rather than excited, like CH+ in collisions with H, H2 
or electrons. We take the effects of excited formation 
and des truc tion into account by the terms Fi and UiDi 
in Eq. IBll Due to the short chemical time-scales, we 
assume the abundance to be constant. At formation, we 
assume the level population to foll ow a Boltzmann dis- 
tribution with temperature Tform (|Stauber fc Brudereii 

^ www.strw.Ieidenuniv.nl/~ewine/photo 
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120091 Ivan der Tak et all [20071) . The destruction is as- 
sumed to be independent of levels. Hence, 



F = Y,F,^Y.Pa, 



, i^ i I ^ -'- ioTTCL 



/Q(7form) 



= ^ ?ij A = D^n,= Dn{GR+) , 



(B3) 



with the partition function Q{T), the statistical weight 
gi of a level with energy Ei. For the calculation of the 
destruction rate D, we take reactions with H, H2 and 
electrons into account. Thus D — fcin(H) + k2n(ii2) + 
k3n(e~), with the reaction rate coefficients fci, ^2 and 
fca. In this approach, the formation temperature Tform 
remains a free parameter. 

Since the level populations at different positions of 
a model are coupled by the radiation field, the solu- 
tion of a radiative transfer problems can be very time- 
consuming, especially in multi-dimensional geometries or 
for molecules with high optical depth like water. We thus 
use the escape probability approach incl uding continuum 
emission and absorption presented by ^Takahas hi et al.1 
(J1983 ) and approximate the ambient radiation by 

(Jij) w (1 — €ij)B{Tcx,ij,'^ij) + {cij — ilij)B{Tdust, Vij) 
+riijB{TcMB,Vtj) , (B4) 

with the line frequency i>ij and the Planck functions 
B[T, v) for the excitation temperature Tex, the tem- 
perature of the dust Tdust and the cosmic microwave 
background Tcmb- The escape probabilities e^ and rjij 



give the probability for a photon to escape line absorp- 
tion {^ij) or to escape both line and dust absorption 
iV y). They are ca l culate d using Eq. 3.11 - 3.13 given 
in iTakahashi et aD ()1983t ). The necessary optical depth 
of dust and lines is obtained from the summation along 
different rays for different frequencies. 

This approach couples different cells and the prob- 
lem is solved iteratively, similar t o one-zone escape 
proba bility codes (e.g. RADEX, Ivan der Tak et aLl 
l2007t ). To improve convergence, we use an ALI- 
like a cc eleration mechanism (e.g. J Rybicki fc Humniieil 
11991! or iHogerheiide fc van der TagHoOOD . To test the 
code, we have run seve ral benchmark problems, e.g. 
Ivan der Tak et aLl (|2005l ). with results typically within 
30 % to the exact soluti on. A similar method has also 
been successfully used by iPoelman fc Spaan^ (|2005l ) . 

The velocity field and the intrinsic line width enter the 
excitation calculation by the ambient molecular radia- 
tion field. These parameters are however only important 
for optically thick lines and the line shape. Optically 
thin lin es are only l i ttle a ffected. The predicted velocity 
field by iDotv et al.l (|2006l ) (after iMcLaughlin fc Pudritd 
Il997f) for young chemical ages of less than a few times 
10^ yrs is approximately static outside a few 1000 AU 
and we thus assume a static situation. The line width 
mainly due to microturbulen ce is assumed to be 1.6 km 
s^^ (Doppler parameter, e.g. iStauber et al.ll200"7l) except 
in the outfiow wall, where a larger line width of 4.2 km 
s~^ is assumed as indicated by interferometric observa- 
tions (|Bruderer et al.ll2009b( ). 
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